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ABSTRACT 

We analyze samples of 5piteer Infrared Spectrograph (IRS) spectra of T Tauri stars in the Ophiuchus, 
Taurus, and Chamaeleon I star-forming regions, whose median ages lie in the < 1 to 2 Myr range. 
The median mid-infrared spectra of objects in these three regions are similar in shape, suggesting, on 
average, similar disk structures. When normalized to the same stellar luminosity, the medians follow 
each other closely, implying comparable mid-infrared excess emission from the circumstellar disks. We 
use the spectral index between 13 and 31 /im and the equivalent width of the 10 ^m silicate emission 
feature to identify objects whose disk configuration departs from that of a continuous, optically thick 
accretion disk. Transitional disks, whose steep 13-31 ^m spectral slope and near-IR flux deficit reveal 
inner disk clearing, occur with about the same frequency of a few percent in all three regions. Objects 
with unusually large 10 /im equivalent widths are more common (20-30%); they could reveal the 
presence of disk gaps filled with optically thin dust. Based on their medians and fraction of evolved 
disks, T Tauri stars in Taurus and Chamaeleon I are very alike. Disk evolution sets in early, since 
already the youngest region, the Ophiuchus core (L1688), has more settled disks with larger grains. 
Our results indicate that protoplanetary disks show clear signs of dust evolution at an age of a few 
Myr, even as early as ~ 1 Myr, but age is not the only factor determining the degree of evolution 
during the first few million years of a disk's lifetime. 

Subject headings: circumstellar matter — stars: formation — stars: pre-main sequence — infrared: 
stars 



1. INTRODUCTION 

Protoplanetary disks around pre-main-sequence stars 
are known to dissipate over the course of a 
few million years; at an age of ~ 10 Myr, 
very f ew objects are still surrounded by primordial 



disks (iJavawardhana et aP 119991: iHaisch et"al 



Sicilia-Aguilar et al.l 120061 IJavawardhana et al 
Hernandez et all l2006l ). During the first 1-2 Myr, disks 
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already experience significant evolution: dust grains have 
been show n to grow and settle towards the disk mid- 
plane (e.g.. iMiv akc fc Nakagawa 1995; A pai et al1l2005t 
iD'Alessio et aT.! .2006; Eurlan et al.. ,200a l200dr^to be 
partl y processed into crysta l hne f o rm (e.g.. iMeeus et al.l 
20031 iKessler-Silacci et al.l 120061: ISargent et all l2006l 



20091 ). and to be removed from certain radial regions in 
a disk, likely due to the gravitational interactions with 
a companion that fo rmed in the disk, such as a newly 
formed planet fe.g.. ISkrutskie et all 119901 : ICalvet et al.l 
I2005t lEspaillat et al.l l2008l). Bv understanding the pro- 
cesses occurring in a protoplanetary disk, we gain in- 

^ NASA Astrobiology Institute, and Department of Physics and 
Astronomy, UCLA, 430 Portola Plaza, Los Angeles, CA 90095 

2 current address: JPL, Caltech, Mail Stop 264-767, 4800 Oak 
Grove Drive, Pasadena, CA 91109; Elise.Furlan@jpl.nasa.gov 

^ Department of Physics and Astronomy, University of 
Rochester, Rochester, NY 14627; dmw@pas.rochester.edu, 
manoj@pas.rochester.edu, khkim@pas.rochester.edu, bsar- 
gent@pas . rochester . edu, forrest@pas . rochester.edu 

* Department of Astronomy, The University of Michigan, 
500 Church St., 830 Dennison Bldg., Ann Arbor, MI 48109; 
melisma@umich.edu, ccespa@umich.edu, ncalvet@umich.edu, 
lhartm@umich.edu 

Centro de Radioastronomfa y Astroffsica, Universidad Na- 
tional Autonoma de Mexico, Apartado Postal 3-72 (Xangari), 
58089 Morelia, Michoacan, Mexico; p.dalessio@astrosmo.unam.mx 



sights into disk clearing and planet formation mecha- 
nisms; observing such disks at the onset of these pro- 
cesses allows us to discriminate better the degree to 
which each process contributes to disk evolution. 

Close to the Sun at a distance of about 150 pc 
(jBertout et al.l 11999 ^. he three star- forming regions 
whose ages are thought to fall into the 1-2 Myr range: 
Taurus-Auriga, Ophiuchus, and Chamaeleon I. Ophi- 
uchus is the youngest (< 1 Myr; Luhman fc Riek3 
Il999( l and most embedded of the three; its core re- 
gion (L16 88; ^ 300 has a very hig h stellar conc entra- 
tion fe.g. jGreenefc You ng 1992; Bo ntemps et al.| [2001). 
Both Chamaeleon I and Taurus are less extinguished 
and crowded than O phiuchus. Chamaeleon I has a me- 
dian age of 2 Myr (|Luhmanl l2004aD . but a l arge age 
sprea d ranging from < 1 Myr up to 6 Myr (|Luhmanl 
|2007| ). and it is somewhat older than Taurus, whose age 
is about 1 Myr, with a distribution of ages between 1 
and several Myr (iKenvon fc Hartmann|[T995: HartmannI 
[200H iL^A^an et al.ll2003i r The region surrounding the 
compact Ophiuchus core contains several young stellar 
objects (YSOs) that are less obscured and older than 
the Ophiuchus core region (Ichikaw a fc Nishidal Il989l : 
iChen et al.l[l995l : lWilking"er a l. 2005). Since the location 
in t he H-R diagram of objects from this off-core region 
(see lWilking et al]|2005D is simila r to that of members of 
Chamaeleon I fsee lLuhma"nll2007[ ) , these two regions have 
comparable ages (~ 2 Myr), independent of evolutionary 
models. 

The majority of objects in the Ophiuchus, Taurus, and 
Chamaeleon I star-forming regions are T Tauri stars, 
which are initially surrounded by optically thick circum- 
stellar disks that are accreting onto the star; during this 
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stage, they are referred to as classical T Tauri stars. The 
dust in the disk causes an excess at infrared wavelengths, 
which decreases over time as the disk material is dissi- 
pated. The slope n of the spectral energy distribution 
(SED), XFx oc A", measured between 2 and 25 /im, is 
used to classify systems at different evolutionary stages: 
Class II objects have — 2 < n < and are surrounded by 
protoplanetary disks, while Class III objects have —3 < 
n < —2 and thus have little or no circumstellar material 
left (iLadall 19871: lAdams et aLlllQSTl : lAndre fc Montmerld 
|1994| ). Pre-main-seauence stars that do not display any 
accretion signatures are also known as weak-lined T 
Tauri stars. 

The infrared excess from protoplanetary disks is emit- 
ted as optically thick, continuum emission from the dusty 
disk interior, and optically thin emission from dust in 
the disk surface layer, whose prominent signatures are 
silicate emission features at 10 and 18 /zm. Thus, mid- 
infrared spectra are ideally suited to study these YSOs; 
they allow us to assess the degree of dust crystallization, 
grain growth and settling, and disk clearing, processes 
which, at an age of 1-2 Myr, are likely in their initial 
stages. Studying protoplanetary disks in the four re- 
gions mentioned above will lead to a better understand- 
ing of how the circumstellar material evolves and eventu- 
ally dissipates, possibly giving rise to planetary systems 
like our own. 

As part of a large survey of star-forming regions 
within 500 pc, our team has obtained several hundred 
5-36 /im spectra of YSOs in the Taurus, Chamaclcon 
I, Ophiuchus core, and Ophiuchus off-core regions us- 
ing the Infrared Spectrograph^ fIRS: lHouck et al.ll2004D 
on b oard the Spitzer Space Telescope (jWerner et al.l 
20041) . These spe ctra are classified and presented in 
Furlan et afl (|2006f) (Taurus), Manoj et al. ( m prepara- 
tion) (Chamaeleon I), and McClure et al. (in prepara- 
tion) (Ophiuchus). Here we analyze some two hundred 
Class II objects identified in those papers. In §[2]we intro- 
duce our sample, observations, and data reduction and 
processing, in §|3]we compare the median IRS spectra of 
our objects in Taurus, Chamaeleon I, and the Ophiuchus 
core region, in §|4]we analyze indicators of disk structure 
and infer the degree of disk evolution in all four regions; 
finally, we discuss our results in § [5] and give our conclu- 
sions in § ini 

2. OBSERVATIONS 

2.1. Sample Selection and Data Reduction 

Our Taurus, Chamaeleon I, and Ophiuchus targets 
were observed as part of an IRS guaranteed-time pro- 
gram during IRS campaigns 3, 4, 12, 19, 20, 21, 22, 
23, 29, and 30, which were scheduled from February 
2004 to April 2006. The objects^ were originally se- 
lected from infrared catalogs of YSOs in these regions 
that were compiled before the launch of the Spitzer Space 
Telescope in 2003 (mainly the IRAS Faint Source and 
Point Source Catalogs; the DENIS Database, 2nd release; 
the 2MASS Catalog; iWilking et aDfigsl ICreene et all 

The IRS was a collaborative venture between Cornell Univer- 
sity and Ball Aerospace Corporation funded by NASA through the 
Jet Propulsion Laboratory and the Ames Research Center. 

''' A table, which lists object names and coordinates, is available 
in the electronic version of the journal; the first few lines of this 
table are shown in the Appendix. 



1994 iKenvon k Hartrna^ 119951: [Barsonv et all 119971 : 
Cambresv et al.l 119981 : iPersi et al.l l2000( ). The targets 



in the Ophiuchus off-co re r egion were chose n from 
llchikawa fc Nishidal (|1989f ) and lChen et"aI1 (|1995f) : about 
half are located in or close to the L1689 cloud, one ob- 
ject (IRS 60) lies in L1709, and the remaining objects 
are spread over a wider area. With the availability of the 
IRS spectra, several targets that previously had an only 
tentative classification could b e identified as Class II ob- 
jects^ fsee lFurlan et al.l ()2006[ ): McClure et al. (in prepa- 
ration); Manoj et al. (in preparation)). In Chamaeleon 
I, we also included a few very low- mass star s and b rown 
dwarfs (spectra l type > M5) from iLuhmanI (|2004a[ ) and 
iLuhmanI (|2007D . which were observed in IRS campaigns 
32, 33, 41, and 42. Our sample consists of 85 Class II 
objects in Taurus, 69 in Chamaeleon I, 63 in L1688, and 
15 in the Ophiuchus off-core regions. 

All targets were observed with either the two low- 
resolution IRS modules (Short-Low [SL] and Long-Low 
[LL], 5.2-14 fim and 14-38 ^m, respectively, A/AA ^ 
90) or the SL module and the two high-resolution mod- 
ules (Short-High [SH] and Long-High [LH], 10-19 ^im and 
19-37 /Ltm, respectively, A/AA ~ 600). In this way we ob- 
tained the full mid-infrared spectrum from 5 to 40 fim for 
each target. While in Taurus most objects were observed 
in mapping mode with 2x3-step maps on the target, 
in Ophiuchus and Chamaeleon I most objects were ob- 
served in staring mode. The spectra were extracted from 
the Spitzer Science Center's basic calibrated data (BCD) 
products, either pipeline versions S13.2 or S14.0 (S16.1 
for the few objects observed in campaign s 41 and 42), 
and u sing the IDL-based SMART package (jHigdon et al.l 
|2004[ ). More details on our observing modes, as well as a 
description of our d ata reduction steps, can be found in 
iFurl an et al.l ()2006[ ): notes on the reduction of individual 
objects are mentioned in McClure et al. (in preparation) 
for Ophiuchus and in Manoj et al. (in preparation) for 
Chamaeleon I. 

2.2. Extinction Corrections 

Basic target properties, such as spectral type, extinc- 
tion, and multiplicity information, were taken mainly 
from the literature (Tables [T][4|). In some cases, opti- 
cal extinction Ay values were not available or resulted 
in dereddened optical photometry inconsistent with the 
expected photospheric emission. For objects with such 
unknown or uncertain extinction that are located in Tau- 
rus or Chamaeleon I, we computed the extinction from 
the observed V-I, I- J, or J-H co lors by applying Mathis's 
extinction curve (|Mathisl[l990h for an Ry value of 3.1. 
For the former two colors, we adopted in t rinsic photo- 
spheric colors from iKenvon fc HartmannI (|1995l ). while 
for the latter one, we assumed an int rinsic J-H color typ - 
ical for a classical T Tauri star from lMeyer et all (|1997l ) 
(i.e., a typical J-H color defined by the CTTS locus). In 
a few cases we slightly adjusted the derived Ay values to 
obtain agreement between observed and expected pho- 
tospheric colors. For objects located in the more deeply 

* Note that in Ophiuchus, where extinction is typically high, our 
Class II sample consists of systems that were identified based on the 
dominance of disk emission and not necessarily on their observed 
2-25 ^m spectral slope (see McClure et al. (in preparation) for 
details). 
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Fig. 1. — Reddening correction for T42, located in Chamaeleon I. Left: The spectrum for A\/=0 is the original IRS spectrum, while 
the other spectra were dereddened using Mathis's extinction curve and successively larger Ay values (from bottom to top); the actual 
extinction of T42 is Ai/=8.18. Right: The IRS spectrum was dereddened assuming Av=8.18 and the three different extinction curves used 
in this paper; the spectrum in the middle is the one we adopted. The inset shows the ratio of the spectra dereddened with the McClure 
curves and the one obtained from dereddening with the Mathis's curve. 



embedded Ophiuchus region, McClure et al. (in prepa- 
ration) derived uniform extinctions and uncertainties for 
the entire sample from observed colors (I-J, J-H, or H-K) 
by applying Mathis's curve for i?y=5.0 (this larger Ry 
value is mor e appropriate for dense, molecular clouds; 
lMathisiri990f ). These authors also obtained near-infrared 
spectra of most of the off-core targets and of several of 
the more embedded core objects to determine spectral 
types, extinctions, and multiplicity. 

In order to create a consistent data set, we decided 
to adopt the same dereddening prescription for all our 
targets, even though extinction curves likely vary among 
star-forming regions. For small Ay values, extinction 
curves for molecular clouds are similar to those for the 
diffuse interstellar med ium (ISM), bii t for Ay > 12 there 
are c lear deviations (jWhittet et al.l Il988t IChiar et"aLl 
|2007| ) . Applying an ISM extinction curve like the Mathis 
curve to objects suffering from larger extinctions causes 
an overcorrection of the silicate emission feature at 9.8 
/im, resulting in a feature that is sharply peaked at 
that wavelength, and a 20/10 /im band ratio that is 
sma ller than typical value s for amorphous silicates (> 
0.4: iDorschner et al.lll995[ ). For example. Figure [T] (left 
panel) shows the silicate feature of T42, a Class II ob- 
ject in Chamaeleon I. The original IRS spectrum displays 
a depression between 9 and 10 fim caused by extinc- 
tion along the line of sight; when correcting the spec- 
trum for reddening using Mathis's curve, the absorption 
disappears. Using the extinction value determined for 
T42, Ay— 8.2, and dereddening the spectrum, the sili- 
cate emission feature appears more square; applying an 
even larger reddening correction will cause a larger in- 
crease around 9.7 /xm (where the correction applied by 
the Mathis curve is largest) , resulting in a sharper silicate 
feature. 

Studying ba ckground objects in molecular clouds, 
iMcClur j (|200i^) found that already an Ay value of 3 
resulted in a mid-infrared extinction curve that was 
different from the Mathis curve; at Ay valu es larger 
than a bout 9, the extinction curves converg ed. iMcClurd 
((2009') derived two new extinction curves for the Ay > 3 
regime, one for 3 < Ay < 9 (0.3 < Ak < 1) and one for 
Ay > 9 (the Ay values quoted here assume Ry— 3.1). 
The effect of applying different extinction curves to the 



IRS spectrum is shown in Figure [T] (right panel); the Mc- 
Clure curves cause a larger correction than Mathis's law, 
with larger differenc es at wave l ength s where ice absorp- 
tions play a role (see iMcClurd ()2009^ for details). 

We therefore adopted three different extinction curves 
to deredden our data, which were applied depending 
on the Ay value of our targets: we used a spline fit 
to the Mathis i?y=5.0 curve for objects with Ay < 3, 
and the two McClure curves for objects with Ay > 3. 
Since the Mathis and McClure laws are given in A\/Aj 
and A\/Ak, respectively, we converted the Ay values in 
Tables [T] to S] to Aj and Ak by dividing by 3.55 and 
9.0, respectively. These three curves were also used to 
deredden fluxes when constructing SEDs with ground- 
based optical photometry from the literature (§ 14. 2|) . 
We note that for IRS spectra the correction for red- 
dening using Mathis's curve for i?y=3.1 or 5.0 is the 
same, since the shape of this extinction curve does not 
vary with Ry ioT n ear- and mid-infrared wavelengths 
(Car delli et al.|[l989l ). Furthermore, most objects in Tau- 
rus and Chamaeleon I suffer from little extinction (the 
median Ay values for our samples are 1.5 and 2.7, re- 
spectively), so the details of the extinction curve have 
little impact on the resulting, dereddened spectra and 
SEDs. On the other hand, the median Ay value for the 
Ophiuchus core is 13.5, and thus uncertainties in the ex- 
act shape of the extinction curve make the dereddened 
Ophiuchus spectra somewhat more uncertain. 

2.3. Sample Characteristics 

Our Taurus and Chamaeleon I samples are repre- 
sentative for the Class II population of their respec- 
tive star-forming regions. From a comparison with 
extensive membership lists (iKenvon &: Hartmannlll995t 
iBriceho et all 120021 : lLuhmanll2004af l. we conclude that, 
for spectral types earlier than MO, our samples in these 
two regions are close to complete, while for later spec- 
tral types (up to M6) we sample at least a few Class II 
objects in each spectral type bin. 

Our sample in Ophiuchus is fairly representative, de- 
spite the high extinction in this region (which renders 
the determination of spectral types difficult); the posi- 
tion of our targets in L1688 is coincident with that o f 
the CO gas that defines the L1688 cloud ()Lorenlll989f) . 
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Fig. 2. — Histogram for the distribution of spectral types among 
the objects entering the median in the Ophiuchus core, Taurus, 
and Chamaeleon I star-forming regions. 

While we do not sample the most deeply embedded ob- 
jects in the Ophiuchus core, the median Ay of 13.5 for 
our L1688 targets suggests that several, if not all, of 
them do belong to the young star-forming core. On the 
other hand, the targets in our Ophiuchus off-core sam- 
ple (especially those not associated with any cloud) are 
likely members of a distinct, presumably older, popula- 
tion in t his region, which has recently been surveyed with 
Spitzer (jPadgett et al.l 120081 ). Our very small sample in 
the off-core areas implies that our targets there are likely 
just the brightest members with an infrared excess. 

3. COMPARISON OF MEDIAN IRS SPECTRA 

In order to generally assess the state of evolution of 
circumstellar disks in our three star-forming regions, we 
computed the median IRS spectrum for the Class II 
objects in these regions. To avoid including objects 
with a large range in luminosities, we restricted the me- 
dian to ob jects with spectral types from K5 to M2 (see 
lD'Alessi~e t al. 199?|). For Taurus, we use d the spectra 
of 55 Class II objects, as in iFurlan et aTl ()2006l ). while 
for Chamaeleon I and L1688, the spectra of 28 Class II 
objects in each region were included. 

The distribution of spectral types of our median sam- 
ples is shown in Figure [21 Despite a median spectral 
type of MO for aU three regions, in L1688 43±12%9of 
objects have spectral types earlier than MO, while the 
this fraction amounts to 31±7% and 29±10% in Taurus 
and Chamaeleon I, respectively. Thus, the median of our 
Ophiuchus core targets is more representative for disks 
around slightly more massive stars than in Taurus and 
Chamaeleon I. 

Before calculating the median for each region, we mul- 
tiplied the dereddened IRS spectrum of each object by 
a scale factor to match the object's dereddened i/-band 
flux to the median H-hand flux value of t he region (near- 
IR ma gnitudes were taken from 2MASS [Skrutskie et al.l 
()2006D ]). The median IRS spectrum was then com- 
puted at each wavelength by using the scaled IRS spec- 
tra. The normalization at H results in an approximate 
normalization to the stellar luminosities, at the same 
time minimizing the effect of uncertain extinctions (see 
iD'Alessio et~al1l999D . The L1688 median carries a larger 

^ Note that the uncertainties in the number fractions quoted 
here and in the following sections are based on Poisson statistics 
and thus represent the standard margin of error. 
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Fig. 3. — The median IRS spectra for Taurus, Chamaeleon I and 
the Ophiuchus core; each median spectrum was normalized at a 
certain, dereddened median flux of the respective region; median 
H-band flux (solid lines) and median J-band flux (dotted lines). 
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Fig. 4. — The median IRS spectra for Taurus, the Ophiuchus 
core, and Chamaeleon I, normalized according to their mean rel- 
ative distances (solid lines) and upper and lower distance limits 
(dash-dotted lines). See text for details. 

uncertainty due to the often uncertain spectral types and 
visual extinctions of the Ophiuchus sources. 

In Figure [3] we show the median IRS spectra for the 
Class II objects in Taurus, Chamaeleon I, and the Ophi- 
uchus core region. In order to check our results, we also 
computed the medians by normalizing at the dereddened 
J-band flux, which is more strongly affected by extinc- 
tion corrections, but less by potential near-infrared ex- 
cess from the circumstellar disks. The resulting medians 
are very similar to the ones normalized at H. For the 
remainder of this paper, when we refer to medians, we 
mean the medians that were internally normalized at H. 

The fact that the Chamaeleon I and Ophiuchus medi- 
ans are about 60% lower and higher, respectively, than 
the Taurus median can partly be attributed to differ- 
ences in the distance to the various regions . Chamaeleon 
I is at 160-170 pc (Whitte Fetlill [1997; Bcrtout et all 
19991) , as opposed to 1 40 pc for Taurus f|Bertout et aLl 
19991 [Torres et al.l 12007,) and 120 pc for Ophiuchus 
([Loinard et al.[ [20081 ). Just due to this diversity, we 
would expect fluxes in Chamaeleon I to be a factor of 
0.7 fainter and those in Ophiuchus to be a factor of 1.4 
brighter. This effect can be seen in Figure [H where the 
median IRS spectra were multiplied by scale factors ac- 
cording to their relative distance to Taurus, first by using 
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Fig. 5. — Histograms for the distribution of scaled, dereddened J- 
and //-band fluxes of the objects that entered the median calcula- 
tion for the Taurus, Chamaeleon I, and the Ophiuchus core regions. 
The flux values were scaled according to the mean distance of the 
various regions with respect to Taurus. 

their mean distance (140 pc for Taurus, 120 pc for Ophi- 
uchus, 165 pc for Chamaeleon I), then by considering the 
distance range typically quoted in the literature (130-150 
pc for Taurus, 160-170 pc for Chamaeleon I; 115-125 pc 
for Ophiuc hus; Bcrtout et al. 1999: ' Whittet et al.|[l997t 
[Torres et al . 2007; Loinard et al. 200^! 

After accounting for the larger distance, the 
Chamaeleon I median approaches the flux level of the 
Taurus median, but it is still somewhat lower (on av- 
erage, at about 85% of the Taurus median if the mean 
distances are adopted). The distance-scaled median for 
the Ophiuchus core is still about 20% higher than the 
Taurus median. This could imply different intrinsic J- 
and iJ-band fluxes, with median near-IR fluxes high- 
est in the Ophiuchus core, followed by Taurus and then 
Chamaeleon I. 

Figure [5] shows the distribution of dereddened J- and 
iJ-band fluxes in the three regions, scaled according to 
their mean distance relative to Taurus to allow a direct 
comparison of the near-IR flux values. The median J- 
and i?-band fluxes in Chamaeleon I are about a factor of 
0.8 smaller than those in Taurus (the factor decreases to 
about 0.7-0.75 if the two outliers with the highest near-IR 
fluxes in Chamaeleon I are disregarded); therefore stars 
in Chamaeleon I are intrinsically fainter. On the other 
hand, the median J- and iJ-band fluxes of the Ophi- 
uchus core region are higher than those of Taurus, by 
factors of 1.26 and 1.11, respectively (Figured]). Thus, 
the intrinsic near-IR fluxes in L1688 are higher than in 
Taurus, and the higher flux level of the distance-scaled 
Ophiuchus median is likely caused by this difference in 
near-IR fluxes alone. 

By normalizing the median IRS spectra of Taurus, the 
Ophiuchus core, and Chamaeleon I to a common _ff-band 
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Fig. 6. — The median IRS spectra for Taurus, the Ophiuchus 
core, and Chamaeleon I, normalized at the dereddened //-band 
median flux of Taurus. 
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Fig. 7. — The median IRS spectra for Taurus, Chamaeleon I and 
the Ophiuchus core around the 10 yaa silicate emission feature, 
after subtraction and division by a continuum fit. The dashed blue 
line represents the emission feature of small (sub-/tm) amorphous 
silicates; it is the average continuum-subtracted and -normalized 
profile of LkCa 15 and CM Aur, scaled to the mean 9.8 /tm flux of 
the three medians. 

flux, which amounts to normalization to a common stel- 
lar luminosity, all three medians essentially overlap (Fig- 
ure Therefore, the mid- infrared excess emitted by 
circumstellar disks in these regions is similar and propor- 
tional to the stellar luminosity; this is expected in typical 
T Tauri disks, where irradiation by the central star dom- 
inates the disk heating (e.g., D'Alessio et al. 1999). The 
shapes of the medians are remarkably similar; on a rela- 
tive scale, the median fluxes agree, on average, within a 
few %. However, in the 5-8 [xva. range, the Ophiuchus me- 
dian is about 15% higher than the Taurus median; this 
difference almost doubles in the 13-17 /im region, with a 
maximum deviation of 35%. This is likely the result of 
applyi ng the extinction curve for Ay > 9 from iMcClurd 
()2009| ) to 3 /4 of the targets that entered the Ophiuchus 
median; this extinction curve exhibits the largest devia- 
tion from the Mathis curve and causes a larger flux in- 
crease around 15 /im (see the inset in the right panel of 
Figure [1]) . Despite some uncertainty in the applicable 
extinction curve for high-extinction regions such as the 
Ophiuchus core, the small differences we see among the 
medians of the three regions are likely real. 

To compare the 10 /im emission features, we first fit 
a continuum, defined as a 3rd, 4th, or 5th order poly- 
nomial, to the 5.6-7.9, 13.0-14.0, 14.5-15.5, 28.0-30.0 /im 
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wavelength regions of each spectrum with an associated 
spectral type of K5-M2. Then, this continuum fit was 
subtracted from the spectrum, and the result divided 
by the fit, which results in continuum-subtracted and 
-normalized 10 /im emission features, i.e., we derived 
the emission of the optically thin dust alone. The me- 
dian of these 10 /xm features for each region is shown 
in Figure [T] Also displayed in this figure is the aver- 
age continuum-subtracted and -normalized silicate pro- 
file of LkCa 15 and GM Aur, two of the Taurus objects 
with t he most pristine, o r interstellar-like, silicate fea- 
tures (|Sargent et al.ll2009f ). 

From Figure [7l it is clear that all three median 10 
fim features are wider than expected for small (sub-/im) 
amorphous silicate grains, as are present in the inter- 
stellar medium, suggesting that grain growth occurred 
in most disks of all three regions. The medians of both 
L1688 and Chamaeleon I display a 10 fiui feature that is 
somewhat more structured than the Taurus median and 
also wider on the short-wavelength side. This could be 
an indication of crystalline silicates, which have narrower 
and more co mplex emission features than amorphou s 
silicates (e.g., iFabian et all 120011 : ISargent et all l2006f ): 
in particular, enstatite and silica have strong features 
around 9 ^m pager et al.lll998t IWenrich fc ChristensenI 
[1996) . The fact that the long- wavelength wing of the 10 
/im feature is not wider for L1688 and Chamaeleon I as 
compared to Taurus suggests that objects in these two 
regions do not have, on average, larger grains than their 
counterparts in Taurus. These results are just tentative, 
since only a model of the silicate profiles will reveal the 
fractions of large and of crystalline grains, but we be- 
lieve that the differences we see in the silicate profiles 
are real, given that the extinction curves we applied to 
our spectra do not introduce artifacts at 10 fj,m. 

Figure [8] shows the median spectra scaled to a com- 
mon H-hand flux (the median of Taurus) , together with 
the quartiles, which delineate the range where 50% of 
spectra lie. For all three regions, the quartiles span a 
comparable range in flux values, and medians and lower 
quartiles follow each other closely over nearly all wave- 
lengths. In the Ophiuchus core region, the lower quartile 
almost matches that of Taurus, while the upper quar- 
tile is generally a bit higher in flux and decreases more 
sharply beyond about 22 /im. For Chamaeleon I, both 
the lower and upper quartiles agree very well with the 
Taurus ones over the 9-22 /im wavelength region. Be- 
yond 22 /xm, the quartiles of the Ophiuchus core and 
Chamaeleon I define a somewhat smaller range in flux 
values around the median than the Taurus quartiles. In 
Taurus, the lower quartile decrease more steeply than the 
upper quartile beyond about 25 fim. Overall, the lower 
quartiles are closer to the medians than the upper quar- 
tiles, signifying that fainter disks span a smaller range in 
infrared excess than the brighter disks. 

4. COMPARISON OF SPECTRAL INDICES AND SILICATE 
FEATURE STRENGTH 

4.1. Definitions and Models 

Besides calculating the median mid-infrared spectra, 
which gauge the typical distribution of dust in circum- 
stellar disks, we computed two quantities that probe dust 
evolution: the spectral index between 13 and 31 /im, 
'^is-aij and the equivalent width of the 10 fim silicate 
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Fig. 8. — The median IRS spectra (solid lines), normalized at 
the dereddened H-hand median flux of Taurus, and their quartiles 
(dashed lines) for Taurus, Chamaeleon I and the Ophichus core. 



emission feature, EW{^O^lm) fsee lWatson et al.ll2009|) . 
The former quantity, formally expressed as 



^113-31 



log(A3lfA3i) -log(Ai3JAi3) 
log(A3i) - log(Ai3) 



(1) 



measures the slope of the spectral energy distribution 
(SED), i.e. XFx oc A". We chose 13 and 31 /im as the 
endpoints for our slope determination, since these wave- 
lengths are dominated by the longer-wavelength contin- 
uum emission from the optically thick disk. The fluxes 
at 13 and 31 /im were measured by averaging the flux 
in the 12.8-14.0 /im and 30.3-32.0 /im region, respec- 
tively; these bands exclude the broad silicate emission 
features at 10 and 20 /xm gen erated in the op t ically thin 
disk atmosphere. As show n in lD'Alessio et al.l (|2006[ ) and 
iFurlan et all (|2005l . l2006f ). a steeper SED is indicative of 
increased dust settling, since a less flared disk intercepts 
less radiation from the central star, resulting in decreased 
heating and thus fainter continuum emission. 

In order to measure the equivalent width of the 10 /ini 
feature for each object, we flrst defined the continuum 
as a polynomial anchored in the 5.6-7.9, 13.0-14.0, 14.5- 
15.5, and 28.0-30.0 /im wavelength regions. Depending 
on the object, we adopted a third-, fourth-, or fifth-order 
polynomial to represent the continuum as a smooth curve 
extending from the wavelength region below 8 /im to 
that beyond 13 /im. Then we subtracted this continuum 
from the observed spectrum and divided the result by the 
continuum; this normalized, continuum-subtracted spec- 
trum was finally integrated from 8 to 13 /im to yield the 
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Fig. 9. — Spectral index between 13 and 31 fim versus the equivalent width of the 10 /xm silicate emission feature for accretion disk 
models with mass accretion rates of 10"^, 10~*, 10~^, and 10"^'' Mq yr~^ {panels from top left to bottom right, respectively). The degree 
of dust settling increases according to symbol size. The colors code the inclination angle of the disk, as indicated in the label of the first 
panel. 



equivalent width of the 10 /im feature: 



EWilOfim)= I E^_I^i^d\ (2) 

-^A. cont 



This definition is actually the exact negative of the stan- 
dard definition of equivalent width, which is positive for 
absorption lines and negative for emission lines. How- 
ever, in order to facilitate representations of the equiv- 
alent width, we chose the above definition. The shape 
of the adopted continuum is often the main source of 
uncertainty of i5VF(10/im); we took this into account 
when determining the overall uncertainty of our calcu- 
lated EW{\QiiTa) values. 

The EW{\QiJLm) is a measure for the amount of op- 
tically thin dust per area of optically thick disk. In a 
typical accretion disk, the optically thin emission is gen- 
erated by smal l dust grains {S , 5 /im) in the disk at- 
mosphere (e.g.. iD'Alessio et aliri999( ). It is the emission 
from these dust grains, and the energy released by accre- 
tion, that heats the lower disk layers; thus, the optically 
thin emission is tied to the optically thick continuum 
emission from the disk interior, and it will depend on 
the overall structure of the disk. As a disk evolves, dust 
grains are thought to grow and settle towards the disk 
midplane; the disk becomes less flared (smaller ni3_3i), 
and the strength of the silicate feature decreases (e.g., 
DuUcmond & Dominik 2004, 2005; D'Alessio ct al. 2006; 
Kessler-Silacci et al.ll200a I Watson et al.ll200a) . 

We determined typical ranges for ni3_3i and 
EW{lQ^va) by computing these quantities for a grid of 



accretion disk models around a 0.5 A/0 star, following 
the same procedures as for the data. We only slightly 
modified the anchor regions for the continuum (7.0-7.9, 
13.5-14.5, 28.0-30.0 /im) due to the different wavelength 
grid of the models. These models were calculated ac- 
_cording the methods of D'Alcssio ct al. (2006) (see also 
Espaillat (2009)); the vertical disk structure was derived 
self-consistently, resulting in a fiared disk with a ver- 
tical inner rim at the dust sublimation radius (located 
where the dust reaches 1400 K). The dust in the disk 
was assumed to be composed of silicates and graphite 
(jPraine fc Led Il984f ) and to scatter isotropically. The 
degree of dust settling was parameterized with e, which is 
the ratio of the adopted dust-to-gas mass ratio in the up- 
per disk layers relative to the standard dust-to-gas mass 
ratio of 1/100. A small value of e implies depletion of 
small grains in the upper disk layers and an increase of 
larger grains close to the disk midplane. 

The results of our ni3_3i and EW{1Q iim) calculations 
for the grid of accretion disk models are shown in Figure 
[9l For the models with a mass accretion rate of 10""^ to 
10~* Mq yr~^, the spread in ni3_3i is mostly attributed 
to dust settling, with the more settled disks having a 
steeper SED slope. Dust settling also causes most of 
the spread in EW{10 ^m) for the 10~® Mq y's:~^ models, 
while for the higher mass accretion rate, this spread is 
mainly caused by different inclination angles. We note 
that at larger inclination angles (> 75°; not shown), the 
silicate feature turns into absorption, especially for the 
more fiared disks (i.e., e ~ 1), resulting in a small and 
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Fig. 10. — Spectral index between 13 and 31 fim versus the equiv- 
alent width of the 10 fim silicate emission feature for the Ophiuchus 
core region. The polygon specifies the region of typical accretion 
disks as derived from models. A few of the "outliers", which are 
discussed in the text, are labeled. 

eventually negative 10 /im equivalent width. 

We obtain a somewhat smaller spread in EW{10fim) 
for the models with a mass accretion rate of 10~^ to 
10~^° Mq yr~^ (Figure[9]), but the range of ni3_3i values 
is larger. The spread in spectral indices can be attributed 
to both dust settling and varying inclination angles. The 
more settled, highly inclined model disks approach the 
spectral index of -4/3, which is the value for an infinite, 
geometrically thin, optically thick disk. 

The less settled model disks have the largest values 
for ni3„3i and EW{10 /Ltm) when their mass accretion 
rates amount to 10~* to 10~^ Mq yr~^, which are repre- 
sentative figures for sta r-forming regions such as Taurus 
([Hartmann et al.l 11998'). In general, when dust grains 
grow and settle, the EW{10fim) decreases; even though 
the presence of larger grains (> 1-2 /im) in the optically 
thin surface layer would increase the width of the silicate 
emission feature, its amplitude would be reduced such 
that the EW{10iim) still decreases. 

The range of disk models defines a region in the ni3_3i 
versus EW{10fi-m) plot which is occupied by typical ac- 
cretion disks. This region extends from -1.33 to about 
1.2 in ni3_3i and from close to to 2-3 (with higher val- 
ues for more flared disks) in EW {10 /im); it is shown as 
the gray polygon in Figures [10] to [131 These latter fig- 
ures show the spectral index from 13 to 31 /xm and the 
equivalent width of the 10 /im feature calculated for the 
Class II objects in the Ophiuchus core region, Taurus, 
Chamaeleon I, and in the Ophiuchus off-core region, and 
are discussed below. 

4.2. Class II Objects 

As is apparent from Figures [10] to [121 niost disks in 
L1688, Taurus, and Chamaeleon I have ni3_3i < (see 
also Fig. [Tl]l : the fraction is essentially the same in Tau- 
rus (80±10%) and L1688 (79±11%) and very similar in 
Chamaeleon I (71±10%). This indicates that the disk 
structures of typical objects in these three regions are, 
on average, comparable, as was gauged from their me- 
dians (§ [3]), and that sedimentation starts early in the 
lifetime o f a protoplanetary disk, v alidating model pre- 



197a iChiang fc GoldreichI 
2005( 1. On the other hand. 



dictio ns (Goldreich & Ward] 
Il997t [Pullemond & Domini^ 
about half of our small sample of the Ophiuchus off- 
core region, which is similar in age to Chamaeleon I, 
has ni3-3i > (Fig. [TS]). Most of these objects typically 
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Fig. 12. — Same as Figure [TOl but for the Chamaeleon I region. 



DoAr 28 
♦ 



16201-2410 



Ophiuchus 
(off core) 

16126-2235 




16225-2607 



5 6 7 8 9 
EW(10 /xm) [fj.m] 



10 11 12 13 



Fig. 13. 
region. 
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have large 10 /im equivalent widths, too, suggesting an 
unusual disk structure. 

In fact, each of these regions contains a substantial 
number of objects with i?W^(10/im) values larger than 
expected from typical accretion disks (see Fig. [14] for 
the distribution of EW {10 fim) values). When count- 
ing the number of objects outside the region defined by 
the polygon in Figures [TOl to [T3l relative to the total num- 
ber of objects (not including objects whose EW{10fim) 
is negative due to poorly defined continua and silicate 
features), we derive that 35±7%, 27±6%, 29±6%, and 
60±20% of Class II objects in the Ophiuchus core, Tau- 
rus, Chamaeleon I, and Ophiuchus off-core, respectively, 
have unexpectedly large silicate equivalent widths and/or 
larger ni3_3i values. Considering the smaller sample 
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Fig. 14. — Histograms for the distribution of rii3_3i and 
EW{10fj,m) values for Taurus, Chamaeleon I, the Ophiuchus core 
and the Ophiuchus off-core regions. 

size, the Ophiuchus off-core region has a marginaUy 
larger fraction of unusual objects than the other three 
regions. There is a lack of objects with very large ni3_3i 
values in the Ophiuchus core region. We note that some 
of the uncertainties in EW{10iim) are relatively large; 
they are the result of ambiguity in the underlying con- 
tinuum, which generally affects objects with strong 10 
fim emission more, and of uncertain extinctions, but the 
latter to a minor extent and mostly for objects in L1688. 

4.2.1. Outliers m ni3_3i 

We define the objects with steep spectral slopes from 

13 to 31 fxm (ni3_3i > 1) as transitional disks. Despite 
being Class II objects, their SED rises sharply beyond 
about 10 /xm, which is attributed to the presence of an 
inner disk wall bounding an inner disk hole whose deficit 
in small dust grains cau ses a decrease in near-infrared ex- 
cess emission jForrrcst et al.ll200 4': 'D'Ale ssio et al.ll2005l: 
ICalvet et al][2005- .Espaillat et al., ,2007a) . A close, stel- 
lar companion could be responsible for the disk hole, and 
a luminous, embedded companion, whose SED peaks in 
the mid- to far-in frared, could create a transitional disk 
appearance (e.g.. rDuchene et al.ll2003[ ). Therefore, inter- 
pretation of transitional disks in terms of disk evolution 
first requires identification of any close companions. 

Ophiuchus core. Only one object associated with 
L1688 has an unusually large ni3_3i value and reduced 
near-infrared excess emission: SR21 (see Table[T]and Fig- 
ure [T5|l . Its IRS spectrum rises sharply beyond about 

14 /im; in addition, its silicate emission feature is su- 
perposed with PAH emission features, which are typi- 
cal for disks arou nd stars of earlier spectral type (e.g., 
iGeers et all 120061 : SR 21 has a spectral type of F4), and 
therefore the 10 /im feature strength cannot be deter- 
mined accurately. SR21 is a 6.4" binary; the IRS spec- 



trum contains the flux of both components, but the sec- 
ondary is unli kely to contribute significantly at infrared 
wavelengths (jPrato et al.ll2003l ). The near- infrared ex- 
cess seen in this object rules out a fully evacuated in- 
ner disk. Since Ha meas urements indicate t hat SR21 is 
not accreting any more (jMartm et al.lfT998f l. the inner 
disk is likely domi nated by dust, not gas. Models by 
I Brown et al.l (|2007D require a gap between 0.45 and 18 
AU to reproduce the observed SED of SR21, with op- 
tically thick disk material remaining in the inner disk. 
Recent high-resolution, sub-millimeter images of this ob- 
ject reve al an inner cav i ty of ~ 37 AU in the dust con- 
tinuum (j Andrews et al.ir2009t ). confirming the presence 
of the inner disk clearing. 

Taurus. In Taurus, four objects have ni3_3i values 
that are well above those found in typical accretion disks: 
CoKu Tau/4, UX Tau A, GM Aur, and DM Tau (Table 
[2] and Figure [T5|) . All have been previously identified 
as transitional disks with different degrees of inner disk 
clearing, and they have been studied in det ail by applying 
SED models (D 'Ales sio et a l. 2005; Calvet et al.l l2005l : 
lEspaillat et all l2007bl) . The disks of CoKu Tau/4 and 
DM Tau are truncated at 10 and 3 AU, r espectively, and 
the inner regions are de pleted in dust (jD'Alessio et al.l 
[20051: ICalvet et al.][2005[ ). GM Aur and UX Tau A have 
inner disk holes at 24 and 56 AU, respectively, but opti- 
cally thin dust within a few AU (GM Aur) or an op- 
tically thick ring of material at the dust sublimation 
radius (UX Tau A) (jCalvet et al.l l2005t lEspaillat etlll 
l2007bl ). Recent submillimeter and millimeter interferom- 
eter maps of GM Aur co nfirm the inner disk hole of ~ 20 
AU ( Hughes et aLll2009f) . Except for CoKu Tau/4, these 
objects are accreting and do not have close companions 
with mass rati os larger than 0.1 oyer a 20-160 mas sep- 
aratio n range teen von et "allll998l: IWhite fc Ghej|2001l : 
Ireland_&_Kraus 200l). Re cent work suggests that CoKu 
Tau/4 is a close binary (jlreland fc KrausI l2008f l. and 
therefore its inner disk clearing is not caused by disk 
evolution. This leaves three "real" transitional disks in 
Taurus. 

Chamaeleon I. In Chamaeleon I, there are five objects 
with a deficit of near-infrared excess emission and large 
ni3_3i values: CS Cha, T25, SZ Cha, T35, and T54 (Ta- 
ble |3| and Figure [p see also lKim et al]l2009h . CS Cha 
is a transitional disk similar to GM Aur, with some op- 
tically thin dust in the inner disk region ijEspaillat et al.l 
|20073)- Despite the re cent discovery that it is a close 
binary system (< 0.1"; iGuenther et al.ll2007f). the fact 
that it s inner hole is ~ 43 AU in size (|Espaillat et al.l 
l2007a^ suggests that the binary might not be the only 
agent clearing out the inner disk of CS Cha. T25 and 
T54 somewhat resemble CS Cha; their emission is pho- 
tospheric below about 10 /im and rises beyon d 15 /im. 
T25 is a single star and lik ely not accreting (jLuhmanI 
l2004at iLafreniere et al.ll2008[ ). implying that disk evolu- 
tion might be respo nsible for the inner disk hole. T54 
is also not accreting (jFeigelson fc Krisill98 9'): instead of 
silicate emission, it displays an 11-13 fim PAH complex, 
and its spectral slope beyond 13 /zm is shallower than 
that of the other transitional dis ks. Since it has a com- 
panion at a separation of 0.27" (jChez et al.lll997l ). the 
inner disk is likely cleared by the binary. Both SZ Cha 
and T35 have some n ear-infrared excess emi s sion and 
are also still accreting (jGauvin fc StromI 119921 : iLuhmanI 
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Fig. 15. — The dereddened SEDs of objects with unusually large ni3_3i values (i.e., larger than expected from accretion disk models); 
the photospheres are represented by Kurucz model atmospheres with solar metallicity, log(g)=3.5, and different effective temperatures: 
Tcff=6500 K (SR21), Toff=3625 K (CoKu Tau/4, T25), Tcff=3750 K (DM Tau), Tcff=4750 K (GM Aur), Toff=4375 K (UX Tau A, 
DoAr 28), Toff=4250 K (CS Cha), Toff=5250 K (SZ Cha), T<,ff=3875 K (T35), Toff=5500 K (T54), and Tcff=6000 K (16201-2410). The 
photometry was adopted from the literature, and the data were dereddened as explained in the text. The dashed gray lines represent 
the medians of the Ophiuchus core (for SR21), of Taurus (for CoKu Tau/4, DM Tau, GM Aur, UX Tau A), and of Chamaeleon I for 
the remaining objects, normalized at the H-band flux of each object. The Chamaeleon I median was used also for the Ophiuchus off-core 
targets, since these two regions are of similar age. 



[20043), suggesting the presence of optically thick ma- 
terial in the in ner disk, as is the case for UX Tau A 
(lEspaillat et al. 2007bl. Whil e T35 is probably a sin- 
gle star (jLafreniere et al.ll2008( ). SZ Cha has two "com- 
panio ns" , one at 5.3" and the other at 12.5" ()Ghez et al.l 
Il997( ): however, these t wo objects are not confirmed 
members of Chamaeleon I ()Luhmanll2007[ ). Therefore, we 
recognize four transitional disks in Chamaeleon I whose 
inner clearings are likely caused by disk evolution. 

Ophiuchus off-core. Finally, we identify two objects in 
the Ophiuchus off-core region with a large 13-31 /um spec- 
tral index and low near-infrared excess emission: DoAr 
28 and IRAS 16201-2410 (Table Hand Figure [l5|). The 
former object lies 0.6° north of the L1688 cloud and is not 
associated with any of the fila mentary clouds in Ophi- 
uchus (e.g.. iRatzka et al1l2005f l. while the latter one hes 
about 0.5° to the west of the main L1688 cloud. The SED 
of DoAr 28 is reminiscent of that of transitional disks 
with inner regions devoid of dust, like DM Tau or T25, 



but it might have a slight excess in the 5-8 /im region. 
Since its Ha emission indicates that it is still accreting 
(■Cohen fc Kuhil[l979h . the inner disk, while depleted in 
dust, likely contains gas. It does not have a compan- 
ion (the upper limit for the b rightness ratio relat ive to 
the primary is 0.04 at 0.15"; iRatzka et all 120051 ) . and 
therefore its disk structure can probably be seen as an 
indication of disk evolution. The IRS spectrum of IRAS 
16201-2410 shows a strong silicate emission feature in ad- 
dition to the steep rise beyond about 14 /im.. Since, as 
opposed to DoAr 28, it has some excess at near-infrared 
wavelengths, the inner disk is not fully cleared of dust. 
Recent near-infrared observations revealed that it has a 
companion at a projected separation of ~ 2" (McClure 
et al., in preparation); since this companion seems to 
be very embedded, it could be responsible for a (small) 
fraction of the long-wavelength excess. 

Three more objects have large rii3_3i values, but dis- 
play no silicate emission feature: IRS48 in L1688, CU 
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Fig. 16. — The dereddened SEDs of three Herbig Ae/Be stars with large ni3_3i values; the photospheres are represented by Kurucz 
model atmospheres with solar metallicity, log(g)=3.5, and different effective temperatures: T(,ff =9500 K (IRS48), T(,ft=10000 K (CU Cha), 
and Tcff =7250 K (16156-2 358). The photom etry for IRS48 was taken from va rious pub lished catalogs (USNO, 2MASS, c2d, IRAS), that for 
CU Cha was adopted from lHillenbrand et al.l f|1992i ). and for 16156-2358 from lVieira et al., ( ,2003 ); the data were dereddened as explained in 
the text. The dashed gray lines represent the median of Ophiuchus (for IRS48) and of Chamaeleon I (for the other two objects), normalized 
at the H-band flux of each object. 
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Fig. 17. — The dereddened SEDs of objects with very low 7113—31 values (< -4/3); the photospheres are represented by Kurucz model 
atmospheres with solar metallicity, log(g)=3.5, and different effective temperatures: Tcff=4000 K (GY292, F04570-I-2520), Tcff=3625 K 
(SRIO), Tofi=5500 K (SR20), Toff=3875 K (VSSG5), and Tcft=4750 K (T51). The photosphere for Hn 5 is an AMES-Dusty model 
atmosphere with Toff=3300 K and log(g)=4.0. The photometry was adopted from the literature, and the data were dereddened as 
explained in the text. The dashed gray lines represent the medians of the Ophiuchus core (for GY292, SRIO, SR20, VSSG5), of Taurus 
(for F04570-I-2520), and of Chamaeleon I (for Hn5 and T51), normalized at the H-band flux of each object. 



Cha in Chamaeleon I, and IRAS 16156-2358 in the Ophi- 
uchus ofF-core (Figure I16p . Their SEDs show numerous 
PAH features and a steep increase beyond about 15 ^m; 
while CU Cha and IRAS 16156-2358 have significant 
near-infrared excess, IRS48 displays no excess below 2 
fim. These objects are considered Herbig Ae/Be stars, 
given their early spectral type (AO for IRS48, B9.5 for 
CU Cha, and FO for IRAS 16156-2358). 

Models of CU Cha showed that, while the near-infrared 
excess can be attributed to the puffed- up inner rim at the 
dust destruction radius, the rise in the SED in the 20- 
30 ^m region is due to the flaring disk that reemerge s 
from the shadow behind the rim (|Doucet et al.l 120071 ). 
A similar disk structure could apply to IRAS 16156- 



2358; in addition, the steep rise of its SED beyond 15 
/zm could be partly due to a companion 4" to the south, 
which is fainter than the primary at 2.2 /im (2MASS, 
ISkrutskie et al.|[2006f ) , but whose mid-infrared flux is not 
known. Therefore, despite the large ni3_3i values, CU 
Cha and IRAS 16156-2358 likely do not belong to the 
transitional disk category. 

On the other hand, the lack of near-IR excess emission 
and the presence of strong PAH fe atures i n IRS48 are 
consistent with observations of Gee rs et al.l ([2007), who 
resolved a gap with a radius of 30 AU in their 18.7 /im 
image of this object, but found that the PAH emission 
filled the gap region. Therefore, IRS48 seems to have 
a transitional disk structure, with a lack of small dust 
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grains in the inner disk regio ns, but with some hot dust 
still present close to the star (|Geers et alJl2007[ ). 

The frequency of disks in a transitional stage, based 
on rii3_3i values, is comparable in all four regions and 
amounts to a few %. It is most uncertain in the Ophi- 
uchus off-core region with 13.3±9.4%, since the sample 
size is very small. The fraction of transitional disks in 
L1688 amounts to 3.2±2.2%; we note that L1688 lacks 
any disks with substantially cleared inner regions that 
would result in a severe deficit of infrared excess out to 
~ 8 /im. The transitional disk fraction in Chamaclcon I 
(5.8±2.9%) is similar to that in Taurus (3.5±2.0%); the 
two star-forming environments are also alike. 

When examining Figures [TUl to [T^ and Tables [1] to [31 a 
few Class II objects with particularly low ni3„3i values 
are apparent: GY292, SRIO, SR20, and VSSG5 in the 
Ophiuchus core, F04570-f 2520 in Taurus, and Hn 5 and 
T51 in Chamaeleon I (see Figure[T7]for their SEDs). The 
13-31 /xm spectral index of these objects is actually lower 
than the —4/3 value expected for an infinite, geometr i- 
cally thin, optically thick disk (e.g., Adams et al.lll98"7l ). 
A possible cause for the steep SED slope is the outward 
truncation of the disk by the gravitational interac tion 
of a close companion (lArtvmowicz fc Lubowl 11994 see 
also lMcClure et al.l (|2008[) for a model of SR20V GY292, 
SRIO, F q4570-|-2520, and H r i 5 do not ha.ye known com- 
panion s (jRatzka et al.ll2005l : iFilrlan et al.l[2006l : iLuhmai] 
[2004al) . while SR20 and VSSG5 are sub-arcsecond bina- 
ries, with separations of 0.04" -0.07" (|Ghez et al.lll995( l 
and 0.15" (jRat zka et al."2005V respectively, and T51 is 
a ~ 2" binary (Corrcia et al. 2006). 

4.2.2. OutUers m EW{10fj,m) 

While the objects with steep SEDs between 13 and 31 
/im have disks whose inner regions are depleted in small 
dust grains to different degrees, the objects with large 10 
/Ltm equivalent widths, but more typical 77.13-31 values, 
are more difficult to explain. The optically thin emis- 
sion is enhanced above levels found for accretion disks in 
hydrostatic equilibrium. Therefore, the ratio of the pro- 
jected area of the optically thin medium to that of the 
optically thick region must be larger than in a typical 
accretion disk. This can be achieved by either increasing 
the amount or emitting area of the optically thin dust, 
or by decreasing the continuum emission generated by 
the optically thick disk regions, but leaving the optically 
thin emission unchanged. 

The most prominent outliers in terms of EW{10iim) 
are found in Chamaeleon I and the Ophiuchus off-core: 
T56, ISO 91, and B43 in the former, and IRAS 16126- 
2235, 16225-2607, and ROX 44 in the latter (see Figure 
[TS]) . Of the three off-core targets, only ROX 44 is associ- 
ated with a cloud, L16 89. Except for IRA S 16126-2235, 
which is a 1.9" binary () Jensen et al.ll200-S ). these objects 
do not have known companions. Their SEDs display a 
deficit in near-infrared excess emission, a strong silicate 
emission feature at 10 /im, and a rising SED beyond 15 
/im, somewhat reminiscent of transitional disks (Figure 
[20]) . In fact, their ni3_3i values, while lower than those 
of transitional disks, lie in a range rarely occupied by typ- 
ical accretion disks. The decreased levels of near-infrared 
excess emission indicate that their inner disks must be 
depleted in dust, and the rise in the SED at 15 fim might 
suggest the presence of some type of inner disk wall. 
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Fig. 18. — The most prominent outliers in terms of EW{10fj,m) 
in Chamaeleon I and in the Ophiuchus ofT-core region, compared 
to the median of Chamaeleon I. 




5 6 7 8 9 10 15 20 25 30 35 

Wavelength (/j.m) 



Fig. 19. — The most prominent outliers in terms of EW{lOfj,m) 
in Taurus and in the Ophiuchus core region, compared to the me- 
dian of Taurus. 

In the Ophiuchus core region, IRS3, GY195, SR9, IRS 
51, and GY292 are notable outliers in EW{10 fj,m), while 
in Taurus LkCa 15, RY Tau, and MHO 3 share this char- 
acteristic (Figure [111). While GY195, GY292, SR9 and 
LkCa 15 have decreased near- infrared excess emission, 
similar to the outliers in Chamaeleon I and the Ophi- 
uchus off-core, IRS3, IRS51, MHO 3, and RY Tau emit 
above median levels for their respective star-forming re- 
gions over the entire IRS wavelength range (Figure [^Hl)- 

IRS3 is a 0.66" binary with a flux ratio of 0.3 at 2.2 
/im (Ratzka et al.l l2005f l , but it is not known to what 
extent the secondary contributes to the infrared excess. 
SR9 also has a companion at a separation of 0.64", but 
i ts flux ratio of 0.0 6 relative to the primary at 2.2 /im 
(Rat zka et al.l l2005l ) renders it unlikely to significantly 
contribute at infrared wavelengths. IRS51 is a binary, 
too, whose secondary at a di stance of 1-65" is ve ry faint 
(flux ratio of 0.04 at 2.2 /xm; iRatzka et alll2005| ). Inter- 
estingly, GY292 has both a large EW{10fj.m) and a 13- 
31 /im spectral index somewhat below -4/3; therefore, its 
outer disk might be truncated, while its inner disk shares 
some of the properties of the other EW{10 /im) outliers. 

RY Tau is the only object among the sample of 
EW{10 ^m) outliers with an early spectral type (Gl; see 
Table [21). It might have a very close comp anion, at a sep- 
aration of at least 24 mas (iBertout et al.|[r999). Based on 
a dip in its mid-infrared SED. lMarsh &: Mahoiievl (|1992f ) 
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Fig. 20. — The dereddened SEDs of objects with v ery prominent 10 /im emission features; the SED of one additional object that also 
belongs to this category {GY292) is shown in Figure [TTl The photospheres for all objects except B43 are represented by Kurucz model 
atmospheres with solar metallicity, log(g)=3.5, and different effective temperatures: Tcff=3500 K (16126-2235, ISO 91, GY195), Tcff=5250 
K (ROX 44), T^ff =4000 K (16225-2607, MHO 3), Tcfl=3750 K (T56), Tcfj=4375 K (LkCa 15, SR9), T^ff = 6000 K (KY Tau), T^g =3625 
K {IRS3), and T<,ff= 4250 K (IRS51). The photosphere for B43 is an AMES-Dusty model atmosphere with Teff=3500 K and log(g)=3.5. 
The photometry was adopted from the literature, and the data were dereddened as explained in the text. The dashed gray lines represent 
the medians of Chamaeleon I (16126-2235, ROX 44, 16225-2607, T56, ISO 91, B43), of Taurus (LkCa 15, RY Tau, MHO 3), and of the 
Ophiuchus core (IRS3, GY195, SR9, IRS51), and normalized at the H-band flux of each object. 



inferred that a gap could be present in its disk. On 
the other hand, recent models suggest that the SED 
and mid-infrared visibilities of RY Tau can be repro- 
duced with an accretion disk surrounded by an opti- 
cally thin envelope, p ossibly replenished by a disk wind 
(jSchegerer et al.ll2008 f). 

Lk Ca 15, which is likely a single star (|Leinert et al.l 
Il993f l. has been recently modeled as a disk with a ~ 46 
AU inner hole that is not fully cleared out. An opti- 



cally thick inner disk wall at the dust sublimation radius 
and an optically thin inner region extending out to a 
few AU were required to reproduc e the near-infrared ex- 
cess a nd 10 fim silicate emission (jEspaillat et al.|[20Q7bl 
120081 ) . High- resolution interferometer images at millime- 
ter wavelengths resolved th e inner ~ 50 AU h ole of LkCa 
15 in the dust continuum (iPietu et al.l 120061 ). validating 
the SED model. It is possible that the other objects with 
an unusually large EW {10 /im), decreased near- and mid- 
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infrared excess, and steep SED rise beyond ~ 15 /im have 
such a disk gap filled with small, optically thin dust. In 
fact. [Andrews et"al] (|2OO 90 recently imaged an inner hole 
of ~ 33 AU in the disk of ROX 44 using high-resolution 
sub-millimeter observations; they also note that, given 
that this object is still accreting and has a notable in- 
frared excess, material must be present within the cavity 
detected in the millimeter continuum emission. 

The fraction of objects with unusually large 
EW {10 fj,m) , but typical ^13-31 values, amounts to 
30±7%, 21±5%, 19±5%, and 47±18% in the Ophiuchus 
core, Taurus, Chamaeleon I, and Ophiuchus off-core 
region, respectively. Thus, Ophiuchus has a marginally 
larger fraction of objects with an uncommon disk 
structure. Even though the Ophiuchus core suffers from 
larger extinction as compared to the other regions, we 
believe that the use of the McClure extinction curves for 
Ay > 3 does not overcorrect the 10 /zm silicate emission 
feature, and therefore our derived values of EW{10 /im) 
are valid. Also, of the seven Ophiuchus off-core outliers, 
only two have Ay > 9, and the four objects with the 
largest EW {10 fim) values suffer from relatively little 
extinction. 

5. DISCUSSION 

5.1. First Steps of Disk Evolution 

Circumstellar disks in the nearby star-forming regions 
of Ophiuchus, Taurus, and Chamaeleon I display clear 
signs of disk evolution. When considering the spectral 
index between 13 and 31 /im and the equivalent width 
of the 10 fim silicate emission feature, outliers that do 
not fall in the region expected for typical accretion disks 
reveal altered, likely more evolved, disk structures. How- 
ever, most of the objects have ni3_3i and EW{10 //m) 
values that are not unusual (see Figures [TUl to [T5|) . When 
comparing them to model calculations, the fact that most 
disks have a negative spectral index suggests dust deple- 
tion factors of 10^ to 10"^ in the upper disk laye rs, as was 
already found for Taurus (|Furlan et al.ll200a ). Hardly 
any objects occupy the region where accretion disks with 
typical mass accretion rates of 10~^ to 10~^ Mq yr~^ 
and e = 1 lie, even in the young Ophiuchus core, imply- 
ing rapid dust growth and settling already at an age of 
< 1 Myr. 

In previous wor k (e.g. , [B ouwman et alj 

200 It Ivan Boekel et al.l l2003t iMeeus et al.. 2002 
Przvgodda et al.l l2003f) the strength of the 10/xm 
silicate feature was analyzed in a different way: it 
was estimated from the peak-over-continuum ratio in 
the 10 fiui region, and it was compared to the flux 
ratio of the normalized spectrum at 11.3 and 9.8 /xm, 
taken as a measure for the shape of the silicate feature. 
For both Herbig Ae/Be stars and T Tauri stars, a 
corre lation between shape and strength is observe d 
(e.g.. Ivan Boekel et al.ir2003HKessler-Silacci et al.ll2006[ ): 
strong features typically have low F11.3/F9.8 ratios, 
which is interpreted as emission from small, amorphous 
silicate grains, as are found in the interstellar medium. 
Conversely, a large F11.3/F9.8 ratio (and weak silicate 
feature strength) suggests the presence of larger and/or 
more crystalline grains, i.e., more processed dust. 

Therefore, the large 10 /im equivalent widths we ob- 
serve in some of our spectra correspond to large peak- 
over-continuum ratios; they could be explained as objects 



with particularly large mass fractions of small grains in 
their disk surface layers. However, this interpretation 
does not take the optically thick disk component into ac- 
count; if there were more small grains in the disk atmo- 
sphere absorbing stellar radiation, the lower disk layers 
would be heated more, resulting in an increase of con- 
tinuum emission and no substantial change in the 10 /um 
equivalent width. We checked this assumption by calcu- 
lating new accretion disk models, which differ from the 
ones presented in § 14.11 only by the maximum grain size 
in the upper disk layers, 0.05 /im as opposed to 0.25 /im 
in the previous models. For these new models, the 10 
/im equivalent width increased by at most 0.2 /im, prov- 
ing that large EW {10 fim) values cannot be explained 
by a large abundance of very small grains in the disk at- 
mosphere. Also changing the dust composition does not 
significantly increase EW{10 fim); we calculated models 
using opaciti es for amorphous silica tes of olivine com- 
position from iDorschner et al.l (|1995f) . which resulted in 
EW{10fj,m) values of up to 4 /im for the less settled 
models (Espaillat 2 009)- 

As noted bv IWatson et al.l (|2009D . the interpretation 
of peak-over-continuum as different degrees of grain pro- 
cessing neglects the likely diversity in the structure of 
the underlying optically thick disk due to dust settling. 
A more settled disk will have decreased continuum emis- 
sion and typically also less emission from the optically 
thin disk surface layers. Thus, our EW{10fim) values 
and also the peak-over-continuum ratios do not neces- 
sarily probe different grain sizes in the disk atmosphere, 
but rather are a measure for disk structure in terms of the 
ratio of the optically thin and optically thick emission. 

5.2. Formation of Disk Structure 

The spectral index measures the degree of dust set- 
tling in a disk, but is also sensitive to disk inclina- 
tion, especially for low mass accretion rates (< 10~^ 
Mq yr^^). An rii3„3i value larger than expected for 
a disk with no dust settling (i.e., e—1) is indicative of 
a substantial change in the structure of the inner disk; 
all transitional disks studied so far in the Taurus and 
Chamaeleon I regions have large ni3_3i values and inner 
disk clearings ranging f rom a few to several tens of AU 
dP'Alessio ct al. 2005; Calvet et al.ll2005l: lEspaillat et all 
l2007ailb : iKim et al.ll2009f) . 

Close binary companions create inner disk gaps whose 
size depends on the orbital parameters of the binary 
(|Artvmowicz fc Lubowl[T99^ . Such systems would ap- 
pear as transitional disks, but in fact be long-lived cir- 
cumbinary disks. The companion star fraction in the 
young star-forming regions of our study is fairly large, 
^ 50 % (for projected separations in the 15-1800 AU 
range: [Chez et al.lll997l ). but it decreases for smaller sep- 
arations of the binary co mponents (^ 5% for 0.02-0.5 
AU projected separations: lMeloll2003D . More specifically, 
40-50% of the stars in our sample are found in multiple 
systems, but not all objects have been observed at suffi- 
ciently high resolution to determine the presence of com- 
panions within ~ 1". Thus, it is likely that a few of the 
transitional disks are the result of clearing by a central 
binary. In fact, CoKu Tau/4 and CS Cha have recently 
been found to be close binary systems ([Ireland fc KrausI 
I2OO8I: iGuenther et al.l [2007,1 . but it seems that only in 
CoKu Tau/4 are the stars responsible for inwardly trun- 
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eating the disk. 

On the other hand, disks are thought to evolve 
from inside out; processes occur on shorter timescales 
in the inner disk, where den sities are higher and 
orbit al speeds larger (e.g., iDullemond fc DominikI 
120041 [Alexander fc Armitagel l2007t ). Substantial 
dust growth, planet formation, photoevaporation, 
or draining of the inner disk induced by the 
magneto-rotati onal instability (MRI) re s ults in inner 
disk clearing ([Ma rsh fc Maho nevI 119921: ICl arke et al.! 
20011: iD'Alessio et a l. 2005; Alexander fc Armitagc 200l1; 



Chiang fc Murrav-C lav 20071: iKim eF al. 2009). If caused 



by these mechanisms of disk evolution, transitional disks 
are in a temporary state, since processes driving the 
clearing of the inner disk will eventually engulf the whole 
disk. The fact that we see evidence of inner disk clearing 
at the young age of L1688 (< 1 Myr) indicates that disk 
evolution starts early. 

The frequency of transitional disks in the three main 
regions we studied (3-6% of objects at an age of 1-2 Myr) 
implies a timescale of several 10** up to ~ 10^ years for the 
transitional disk phase, assuming disk evolution is con- 
tinuous and that all disks experience this stage in their 
evolution. It is also possible, but unlikely, that only 3-6% 
of disks in a star-forming region appear as transitional, 
with the mechanism that causes the inner disk hole ac- 
tive for a long period (~ 1-2 Myr) of time, such as in 
a circumbinary disk. Disks are known to dissipate over 
time, with clearing processes on the timescale of ^ 10^ 
years, as model calculations for disk clearing by planet 
formation (e.g., Quillen ct al. 2004; Varnicrc ct al. 2006), 
photoevaporation (e.g. Alexander et al. 2006), and MRI- 
activated disk draining ((Chiang fc Murrav-Clav. i2007f) 
have shown. The lifetime of a transitional disk is thus 
shorter than th at of a primordial disk ('~ a few Myr). 
Recent work bv lCurrie fc Kcnyon (2009) suggests a long 
transitional disk timescale, but this is mainly due to the 
different definition of transitional disks; those authors 
consider all "evolved" disks as transitional, while we only 
include disks that show evidence of inner disk clearings. 
By adopting a similar definition of transitional disks as 
we do and using high-quality Spitzer data for Taurus 
and a consistent me thod of estimating disk fractions, 
iLuhman et al.l ()2009D show that transition timescales ap- 
pear to be considerably less than primordial disk life- 
times, substantiating our result. 

The equivalent width of the 10 /im silicate emission 
feature could also suggest changes in the configuration of 
the disk, especially when accompanied by a decrease in 
the near- and mid-infrared emission. A large EW{10 fim) 
implies stronger emission by optically thin dust relative 
to the optically thick disk emission; a gap in the disk, 
filled with optically thin dust, could create this signa- 
ture and would be a precursor of a transitional disk with 
an inner disk hole ()Espaillat et al.ll2007bl . I2008D . On the 
other hand, an object with a dust-free gap would not ex- 
hibit a large EWjlOu m), as is the case for UX Tau A 
(|Espaillat et al.l [20071)1 ). An alternative explanation for 
the EW (10 ^.m) outhers, a s shown by modeling of RY 
Tau (|Schegerer et al.ll2008f ). is the presence of an enve- 
lope of optically thin dust that could be a source for the 
enhanced optically thin emission seen in these objects. 
However, this envelope would have to be replenished in 
small dust grains by some type of flows originating in the 



disk, and it is not known whether sufficiently high mass 
loss rates could be achieved. 

Assuming most objects with large EW{10fim) have 
developed gaps in their disks and that all disks experi- 
ence this evolutionary phase once in their lifetime, the 
larger number of objects presumably with gaps in their 
disks compared to that of objects with inner holes ('^ 
20-30% of objects at an age of 1-2 Myr) implies that 
the stage in which a disk has a gap lasts longer than the 
period of inner disk clearing, from several 10^ up to 10^ 
years. Alternatively, a disk with a gap could survive for 
a shorter period of time, but then experience this phase 
a few times before being dissipated. In addition, even 
if just a few of the disks with large EW{10 /im) in the 
Ophiuchus core were disks with gaps, the presence of 
such disks in this young region would suggest that gap 
formation sets in early. It is interesting that ~ 50% of the 
objects in the Ophiuchus off-core region have unusually 
strong 10 /im emission features; if the gap interpretation 
holds true, it would be a region with a large fraction of 
evolved disks. 

Of the disk clearing mechanisms mentioned above, 
planet formation is the most likely explanation for the 
existence of a gap between inner and outer optically thick 
disk regions. Photoevaporation and MRI-induced inner 
disk draining would not allow for a remnant inn er, op- 
tically thick disk region (see also lEspaillat et al.l [2008) , 
and it is unlikely that dust grain growth would be limited 
to a certain, well-defined radial range in the disk. Since 
a planet opening a gap will likely migrate inwards on 
roughly the viscous diffusion timescale of 10^ initial or- 
bita l periods (e.g. jLin fc Papaloizoulll986l : iNelson et al.l 
I200Q1 . and the duration of the gap phase as estimated 
above could be at least a factor of a few larger, it is pos- 
sible that in some systems we would be witnessing the 
effects of already the second or even later generation of 
planets. This t ype o f scenario has also been proposed by 
IWatson et al.l (|2009D to explain the variety of crystalline 
silicate abundances observed in T Tauri stars in Taurus. 
If we assume that at least some of the objects we studied 
have a disk gap opened by a planet, then these planets 
must have accumulated a mass of at least a few tens of 
Earth masses in < 1 Myr such that they are massiy e 
enough to open up gaps (see, e.g.. iTakeuchi et aT]ll996l ). 
Gravitational instability can create massive planets on 
a timescale of ~ 10^ years (e.g., [Boss 2000), but also 
a migrating protoplanet can grow to a sufhciently large 
mass within ^ 10^ years by accretion of planetesimals 
(jTanaka fc Idal[T999l) . 

A few objects with large ni3_3i values also have large 
10 /im equivalent widths (e.g., DM Tau, CM Aur). The 
large spectral index suggests an inwardly truncated op- 
tically thick disk, with a disk wall at the inner radius. 
For cases where the presence of a near-infrared excess 
indicates that optically thin dust lies in the inner re- 
gions, the unusual strength of the 10 /im f eature can be 
explained by emission from this dust (see ICalvet et al.l 
.2005.) . On the other hand, for objects with no excess 
below about 8 /im and therefore dust-depleted inner re- 
gions, the 10 um emission likely ari s es in the atmospher e 
of the wall (jP'Alessio et al.l 120051 : ICalvet et all |2005l). 
Thus, in a few cases, not only the object's location in 
rii3_3i-£'W^(10/im) space, but also the detailed shape of 
the SED are necessary to interpret disk structure. 
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5.3. Median Infrared Excess 

The interpretation of disk evolution based on the me- 
dian IRS spectrum is compUcated by issues involving 
normalization. If the distribution of source fluxes in 
the near-IR and differences in the distances to the star- 
forming regions are taken into account, then the Class II 
objects in the Ophiuchus core, Taurus, and Chamaeleon 
I have very similar mid-infrared excess emission; the me- 
dians follow each other closely. This indicates that, on 
average, disk structures are similar and that comparable 
amounts of stellar radiation are reprocessed in the disk. 

In particular, the Chamaeleon I star-forming region is 
very similar to Taurus in its disk evolution, as gauged 
from both the median IRS spectra of T Tauri stars and 
the fraction of objects with unusual 13-31 spectral slopes 
and 10 ^m equivalent widths. Agreement between these 
two regions is expected since their molecular cloud envi- 
ronments are comparable and their age difference is min- 
imal, especially cons idering that there is an age spread 
in both regions fe.g. JPalla fc Stable j|2000f) . w hich might 
just b e a result of observational uncertainties (|Hartmannl 
l2Ci0l . 

Based on the distribution of near-IR fluxes, we found 
that the stars in Chamaeleon I are intrinsically fainter, 
while those in L1688 are brighter than their counterparts 
in Taurus. The latter observation can partly be explained 
by the somewhat higher fraction of earlier-type stars in 
our Ophiuchus sample as compared to the other two. 
In addition, since the Ophiuchus core region is younger 
than Taurus, we would ex pect more luminous young stars 
(e.g. jBaraffe et al.ll2002[ ). The age of 0.3-1 Myr for the 
Op hiuchus core was derived assuming a distance of 165 
pc (|Luhman fc Rie ke'lOQi/) ; the ne w distance determi- 
nation of 120 pc by Loinar d et al.l (2008) implies that 
stars are a factor of 1.9 less luminous. Assuming they 
are contracting along Hayashi tracks at roughly constant 
effective temperature {L oc t~^^^: [Hartmannlll998f ). this 
would imply an increase in age by a factor of 2.6. There- 
for e, the median age of ^ 0.3 Myr found for Ophiuchus 
bv lLuhman fc Riekel (|1999( ) would change to ~ 0.8 Myr, 
making it just slightly younger than Taurus. The some- 
what older age of Chamaeleon I, as compared to the other 
regions, could explain the reduced luminosity of its stars. 

The shape of the median 10 /im silicate emission fea- 
ture in the three regions suggests that the disks in L1688 
and Chamaeleon I might have, on average, more crys- 
talline silicates than those in Taurus. Since some pro- 
cessing at high temperatures is required to transform the 
amorphous silicate grains fr om the interstellar me dium 
into crystalline form (e.g., iHarker fc DeschI l2002t l. we 
would expect an older star-forming region, which had 
more time for this processing to take place, to have a 
larger fraction of crystalline grains. On the other hand, 
the larger degree of dust processing in the youngest re- 
gion of our sample, L1688, could be a result of different 
initial conditions in the dense Ophiuchus core region. If 
the mass accretion rates or the masses of the disks were 
higher, disks would likely evolve faster. This could also 
explain the large fraction of settled disks (those with 
n-is-si < 0) in L1688, which is basically the same as 
in the slightly older Taurus region. The fact that our 
Ophiuchus core sample seems to be biased towards some- 
what more massive, young stars supports the idea of 



faster disk evolution, given that there are indications that 
intermediate-mass stars evolve faster than low-mass stars 
(jHernandez et al.ll2005[) . 

6. CONCLUSIONS 

We analyzed the IRS spectra of Class II objects in 
the Ophiuchus core, Ophiuchus off-core, Taurus, and 
Chamaeleon I star- forming regions. Our main conclu- 
sions are as follows: 

• The spectral index between 13 and 31 /im, rii3-3i, 
and the equivalent width of the 10 /im silicate emis- 
sion feature, EW{10^m), are useful indicators for the 
degree of disk evolution. We interpret Class II objects 
with unusually large ni3_3i as transitional disks with 
inner disk clearings, while objects with very substantial 
-^^^(lO fj,m) and decreased near- and mid-infrared excess 
emission could have gaps in th eir disks that are fille d 
with optically thin dust (see also lEspaillat et aT]|2007bD . 
Considering various disk clearing mechanisms, the most 
likely explanation for gaps within optically thick disks is 
planet formation. 

• The Ophiuchus off-core region, which is of similar 
age as Chamaeleon I (~ 2 Myr), has a marginally larger 
fraction of Class II objects whose ni3_3i and EW{10 /im) 
are outside the range found for typical accretion disks, 
as compared to the other three regions. However, the 
sample size is small, and therefore, if these outliers are 
interpreted as objects with an evolved disk structure, we 
can only tentatively conclude that the Ophiuchus off-core 
is the most evolved region. 

• The transitional disk fraction amounts to a few % in 
all four regions; while L1688 lacks any disks with inner 
regions fully depleted in small dust grains, the frequency 
of transitional disks in Chamaeleon I is comparable to 
that of the Taurus region. Assuming disk evolution is 
continuous and affects the structure of the inner disk at 
some point in the lifetime of a disk, we infer that the 
transitional disk phase lasts for several lO"' up to ^ 10^ 
years. 

• The shape and normalized flux levels of the me- 
dian mid-infrared spectrum of disks in Taurus and 
Chamaeleon I are very similar, and the fractions of 
objects with unusually large EW{10fim), but typical 
'^13-31 values, are essentially the same. This sug- 
gests that the similar star-forming environments and age 
ranges of these two regions result in an overall compara- 
ble degree of disk evolution. 

• The median mid-infrared excess of T Tauri stars in 
the Ophiuchus core region is also similar to that in Tau- 
rus and Chamaeleon I. The shape of the median 10 fj,m 
silicate emission feature hints at the presence of larger 
grains and crystalline silicates, and therefore dust grains 
already experienced some processing. This region also 
has a sizable fraction 30%) of Class II objects with 
large 10 /im equivalent widths, which could indicate that 
disk evolution sets in early. 

We are likely witnessing the combined effect of age and 
molecular cloud environment on disk evolution. The in- 
ner disks in the Taurus, Chamaeleon I, and Ophiuchus 
off-core regions, but also in the younger and denser Ophi- 
uchus core, show clear signs of evolution, like the pres- 
ence of larger grains and decreased flaring, indicative of 
dust sedimentation. Some disks also display evidence for 
structure formation, like the development of inner holes 
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TABLE 5 
Target names and coordinates 



Name R.A. (J2000) Dec. (J2000) 



04108- 


h2910 


04 13 57.38 


29 18 19.3 


04187- 


hl927 


04 21 43.24 


19 34 13.3 


04200- 


h2759 


04 23 07.77 


28 05 57.3 


04216- 


h2603 


04 24 44.58 


26 10 14.1 


04303- 


h2240 


04 33 19.07 


22 46 34.2 



Note. — Table 5 is published in its en- 
tirety in the electronic edition of the Astro- 
physical Journal. A portion is shown here for 
guidance regarding its form and content. 



TABLE 6 

Ophiuchus core median (H normalization) 



A 
(^tm) 


(erg/s/cm^) 


(erg/s/cm^) 


(erg/s/cm^) 


5.0 


4.960E-10 


3.300E-10 


7.499E-10 


5.1 


4.561E-10 


2.778E-10 


7.436E-10 


5.2 


4.002E-10 


2.350E-10 


5.905E-10 


5.3 


3.665E-10 


2.123E-10 


5.743E-10 


5.4 


3.509E-10 


2.043E-10 


5.385E-10 


5.5 


3.241E-10 


1.907E-10 


4.998E-10 


Note. 


— Table 6 


is published in 


its entirety in 



the electronic edition of the Astrophysical Journal. 
A portion is shown here for guidance regarding its 
form and content. 



and possibly also gaps. This suggests that grain growth, 
settling, and crystallization occur fast, but the develop- 
ment of radial structure in a disk requires, on average, 
more time and might be a transient phenomenon. Our 
analysis implies that the degree of disk evolution does not 
necessarily correlate with age during the first few million 
years, since it sets in relatively early, and the subsequent 
evolution might be dictated more by the individual sys- 
tem properties than some universal processes happening 
at set time intervals. 



We thank the referee for a thoughtful review that 
led us to improve this paper. This work is based on 
observations made with the Spitzer Space Telescope, 
which is operated by the Jet Propulsion Laboratory 
(JPL), California Institute of Technology (Caltech), un- 
der NASA contract 1407. Support for this work was 
provided by NASA through contract number 1257184 
issued by JPL/Caltech. E.F. was partly supported 



by a NASA Postdoctoral Program Fellowship, admin- 
istered by Oak Ridge Associated Universities through a 
contract with NASA, and partly supported by NASA 
through the Spitzer Space Telescope Fellowship Pro- 
gram, through a contract issued by JPL/Caltech un- 
der a contract with NASA. N. C. and L. H. acknowl- 
edge support from NASA Origins grants NNG05GI26G, 
NNG06GJ32G, and NNX08AH94G. P.D. acknowledges 
grants from CONACyT, Mexico. This publication makes 
use of data products from the Two Micron All Sky Sur- 
vey, which is a joint project of the University of Mas- 
sachusetts and the Infrared Processing and Analysis Cen- 
ter/Caltech, funded by NASA and the NSF. It has also 
made use of the SIMBAD and VizieR databases, oper- 
ated at CDS (Strasbourg, France), NASA's Astrophysics 
Data System Abstract Service, and of the NASA/ IPAC 
Infrared Science Archive operated by JPL, Caltech, un- 
der contract with NASA. 
Facilities: Spitzer(IRS) 



APPENDIX 

Additional tables are provided in the electronic edition. To show their format and content type, the first few lines 
of these tables are printed here. Table 5 contains the target names and coordinates of our sample of Class II objects 
in Taurus, Chamaeleon I, and the Ophiuchus core and off-core regions. Table 6 contains the median IRS spectra and 
their lower and upper quartiles for the Ophiuchus core, Chamaeleon I, and Taurus regions, normalized at the 2MASS 
H-band fluxes. Here, only the first few entries of the Ophiuchus core median and quartiles are shown. These are the 
medians displayed in Figure [3] as solid lines. 
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TABLE 1 

Properties of the targets in the Ophiuchus core region 



Name 


Spectral Type" 


Av'' 


EW{10 tim) 


"13 31 


Median 




(2) 


f3) 


(4) 


(5) 


(6) 


16220-2452 


M3 


3.7 


4.87 


-0.09 




16237-2349 


K5.5 


5.0 


1.67 


-0.17 


X 


B162713-241818 


MO 


13.0 


1.39 


-1.28 


X 


DoAr24 


K5 


4.1 


3.78 


-0.39 


X 


GSS29 


Ml 


11.5 


2.65 


-1.15 


X 


GSS31 


G6 


7.0 


1.58 


-0.50 




GSS37 


MO 


8.7 


1.47 


-1.23 


X 


GSS39 


MO 


17.2 


1.15 


-0.75 


X 


GY144 


M5 


40.3 


3.23 


-0.86 




GY154 


M6 


23.3 


2.63 


-0.27 




GY188 


M3 


40.7 


1.85 


0.15 




GY195 


M3 


24.7 


6.23 


-0.15 




GY204 


M5.5 


3.4 


3.41 


0.22 




GY213 


M4 


27.3 


3.21 


-0.95 




GY224 


M4 


34.2 


3.94 


-0.49 




GY235 


M5 


8.7 


3.02 


0.14 




GY245 




44.2 


3.53 


-0.72 




GY260 


M4 


31.7 


2.31 


-0.18 




GY284 


M3.25 


7.0 


3.60 


0.08 




GY289 


M2 


24.0 


1.54 


-0.20 


X 


GY292 


K7 


13.5 


6.31 


-1.41 


X 


GY301 


K7 


40.4 


3.15 


-0.26 


X 


GY310 


M4 


7.5 


2.92 


0.31 




GY314 


K5 


8.1 


1.08 


-0.20 


X 


GY323 


M5 


32.8 


2.99 


-1.17 




GY326 


M2 


10.2 


2.44 


-0.48 


X 


GY344 


M6 


19.1 


1.35 


0.02 




GY350 


M6 


8.5 


2.53 


-0.94 




GY352 


M5 


20.3 


0.91 


-1.10 




GY371 


M6 


6.3 


0.44 


-0.71 




GY397 


M6 


5.8 


1.18 


-0.76 




HD 147889 


B2 


4.4 


0.02 


-0.11 




IRS2 


K3.5 


7.7 


1.53 


-0.23 




IRS26 


M6 


19.6 


3.41 


-0.78 




IRS3 


M2 


8.5 


6.56 


0.37 


X 


IRS33 


M2 


40.5 


0.99 


-1.08 


X 


IRS34 


MO 


28.7 


0.85 


-0.94 


X 


IRS35 




45.7 


3.65 


-1.03 




IRS42 


K7 


29.0 


3.08 


-0.90 


X 


IRS45 


K6.5 


24.5 


1.59 


-0.87 


X 


IRS47 


M3 


28.0 


1.66 


-1.33 




IRS48 


AO 


15.0 


-0.30 


1.33 




IRS49 


K5 


12.2 


3.67 


-0.11 


X 


IRS51 


K6 


37.6 


5.53 


-0.64 


X 


S2 


MO 


11.8 


0.55 


-1.33 


X 


SR4 


K4.5 


2.6 


2.29 


0.29 




SR9 


K5 


1.9 


5.87 


0.02 


X 


SRIO 


M2 


1.5 


3.03 


-1.34 


X 


SR20 


G7 


7.0 


1.21 


-2.53 




SR21 


F4 


8.0 


4.01 


1.57 




VSS27 


G3.5 


8.1 


1.53 


-0.97 




VSSGl 


MO 


16.5 


0.88 


-1.25 


X 


VSSG5 


MO 


20.3 


1.00 


-1.46 


X 


VSSG25 


M4 


13.6 


3.38 


-0.21 




WLl 


M4 


28.2 


4.51 


-0.19 




WL2 


MO 


34.8 


2.63 


-0.68 


X 


WL3 


M4 


42.0 


1.34 


-0.42 




WL4 


Ml. 5 


21.3 


1.77 


0.06 


X 


WLIO 


K7 


13.3 


1.08 


-0.31 


X 


WLll 


MO 


15.9 


1.97 


-0.41 


X 


WL18 


K6.5 


12.2 


1.40 


-0.34 


X 


WSB37 


M5 


2.9 


1.23 


-0.88 




WSB60 


M4.5 


4.5 


2.01 


0.05 





Note. — Column (1) lists the name of the object, column (2) the adopted spectral type, column (3) the optical extinction Ay 
(assuming Rv=S.l), column (4) the 10 ^m equivalent width {EW {10 j-im)) , column (5) the 13-31 fim spectral index (ni3_3i), 
and the x symbol in column (6) indicates that the object was used in the median calculation. Note that uncertainties in Ay, 
ranging from 1% to 30% of the Ay value (McClure et al., in preparation), account for ~ 25% of the error bars of EW {10 ^im) 
and ni3-3i. 

° Spe ctral types were taken from iBrandner Zinneckeil (|1997f ): iLuhman fc Riekd (|1999f ): IWilking et~all (120051 ): iNatta et all 

(|2006l ): McClure et al. (in preparation). 

^ Av values are from McClure et al. (in preparation). 



Disk Evolution in Three Nearby Star-Forming Regions 



TABLE 2 

Properties of the targets in the Taurus region 



Name 


Spectral Type" 


*■ V 


EW{10 nm) 




Median 




(2) 


(3) 


(4') 


(5) 


(6) 


04108+2910 


MO 


1.4 


0.28 


-0.79 


X 


04187+1927 


MO 


0.0 


0.94 


-0.84 


X 


04200+2759 




0.0 


1.52 


-0.41 




04216+2603 


Ml 


0.0 


1.16 


-0.41 


X 


04303+2240 




11.7 


2.21 


-1.15 




04370+2559 




9.8 


4.06 


-0.51 




04385+2550 


MO 


7.8 


2.48 


0.26 


X 


AA Tau 


K7 


1.8 


1.79 


-0.51 


X 


AB Aur 


AO 


0.3 


4.07 


0.81 




BP Tau 


K7 


1.0 


2.64 


-0.58 


X 


CI Tau 


K7 


2.0 


2.55 


-0.17 


X 


CoKu Tau/3 


Ml 


5.0 


2.66 


-1.19 


X 


CoKu Tau/4 


M1.5 


3.0 


5.37 


2.12 


X 


CW Tau 


K3 


2.8 


1.15 


-0.65 




CX Tau 


MO 


1.3 


2.30 


-0.15 


X 


CY Tau 


K7 


1.7 


0.77 


-0.99 


X 


CZ Tau 


M1.5 


2.4 


2.40 


-1.04 


X 


DD Tau 


M3.5 


1.0 


1.10 


-0.74 




DE Tau 


MO 


1.2 


1.81 


-0.13 


X 


DF Tau 


MO 


1.6 


0.83 


-1.09 


X 


DG Tau 




1.6 


0.13 


0.15 




DH Tau 


MO 


1.7 


3.57 


0.26 


X 


DK Tau 


MO 


1.3 


3.79 


-0.80 


X 


DL Tau 


K7 


1.5 


0.51 


-0.77 


X 


DM Tau 


Ml 


0.7 


5.80 


1.29 


X 


DN Tau 


MO 


0.6 


1.06 


-0.42 


X 


DO Tau 


MO 


2.0 


0.94 


-0.13 


X 


DP Tau 


M0.5 


0.6 


1.70 


-0.33 


X 


DQ Tau 


MO 


1.6 


0.65 


-0.41 


X 


DR Tau 




1.2 


1.23 


-0.40 




DS Tau 


K5 


1.1 


2.35 


-0.96 


X 


F04101+3103 


Al 


1.9 


4.46 


-0.08 




F04147+2822 


M4 


2.5 


2.55 


-1.11 




F04192+2647 




0.0 


1.17 


-0.47 




F04262+2654 




0.0 


0.78 


-0.13 




F04297+2246A 




0.0 


3.93 


-0.41 




F04570+2520 




0.0 


0.94 


-1.60 




FM Tau 


MO 


1.4 


3.13 


-0.16 


X 


FN Tau 


M5 


1.4 


1.84 


-0.05 




FO Tau 


M2 


3.0 


1.40 


-0.32 


X 


FP Tau 


M4 


0.0 


0.94 


-0.10 




FQ Tau 


M2 


1.9 


1.23 


-0.46 


X 


FS Tau 


Ml 


1.4 


1.28 


0.08 


X 


FT Tau 


C 


0.0 


1.86 


-0.46 




FV Tau 


K5 


5.3 


1.47 


-0.66 


X 


FX Tau 


Ml 


2.0 


4.04 


-0.39 


X 


FZ Tau 


MO 


3.7 


1.26 


-0.87 


X 


GG Tau 


MO 


1.0 


2.82 


-0.31 


X 


GH Tau 


M2 


1.0 


1.42 


-0.33 


X 


GI Tau 


K6 


2.3 


2.65 


-0.63 


X 


GK Tau 


MO 


1.1 


4.60 


-0.37 


X 


GM Aur 


K3 


1.2 


5.07 


1.75 




GN Tau 


M2 


3.5 


2.93 


-1.01 


X 


GO Tau 


MO 


2.0 


1.99 


0.03 


X 


Haro 6-13 


MO 


11.9 


3.87 


0.38 


X 


Haro 6-37 


K7 


3.8 


1.38 


-0.85 


X 


HK Tau 


M0.5 


2.7 


1.80 


0.74 


X 


HN Tau 


K5 


1.5 


2.44 


-0.45 


X 


HO Tau 


M0.5 


1.3 


2.46 


-0.63 


X 


HP Tau 


K3 


2.8 


2.43 


-0.01 




HQ Tau 




0.0 


3.05 


-0.64 




IP Tau 


MO 


0.5 


4.18 


-0.09 


X 


IQ Tau 


M0.5 


1.4 


1.75 


-1.00 


X 


IS Tau 


MO 


3.2 


2.30 


-1.21 


X 


IT Tau 


K2 


3.8 


1.12 


-0.87 




LkCa 15 


K5 


1.2 


6.70 


0.62 


X 


MHO 3 


K7 


8.3 


5.93 


0.27 


X 


RW Aur A 


K3 


0.5 


1.25 


-0.54 




RY Tau 


Gl 


2.2 


6.51 


-0.09 




SU Aur 


Gl 


0.9 


4.75 


0.74 




T Tau 


KO 


1.8 


-0.27 


0.55 




UX Tau A 


K5 


0.7 


0.78 


1.83 


X 


UY Aur 


K7 


2.1 


2.05 


-0.04 


X 


UZ Tau 


Ml 


1.0 


2.25 


-0.73 


X 


V410 Anon 13 


M6 


5.8 


1.70 


-0.58 
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TABLE 2 — Continued 



Name 


Spectral Type" 


Ay"' 


EW{IQ ^im) 


»T.13-31 


Median 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


V7in Tail 


Ml 


1.9 


1.43 


-0.70 




V773 Tau 


K3 


2.0 


0.82 


-0.85 




V807 Tau 


K7 


0.6 


0.67 


-0.18 


X 


V836 Tau 


K7 


1.1 


3.34 


-0.45 


X 


V892 Tau 


B9 


8.0 


3.06 


-0.13 




V955 Tau 


K5 


3.7 


1.54 


-0.94 


X 


VY Tau 


MO 


1.4 


2.34 


-0.13 


X 


XZ Tau 


M2 


2.9 


0.66 


-0.32 


X 


ZZ Tau 


M3 


1.4 


1.01 


-0.89 




ZZ Tau IRS 


M4.5 


1.5 


1.58 


0.17 





Note. — S ec the footnote of Table [T] for an explanation of the table content. Uncertainties in Ay, assumed to be 0.3 mag 
(jKenvon fc Har tmann 1995), account for ~ 15% and 1% of the error ba r s of EW(\Oimi) and W13-31 , respectively. 
" Spectral type and extinction information wore taken from Kcnyon ct al. (1990); Hartigan ct al. (1994); Strom & Strom (1994 



iKenvon fc HartmannI (Il995l'l ; iTorres et al ] ( , 1995): Kenyon ct al. (1998); Bric cfio ct al, [199 8); Luhman ( 2 000); White k Ghe3 
I2001 h : iHartiean fc Kenvonl (|2003l ): I White fc BasrinpOOSh : IWhite fc HiUenbrandl (|2004l ): ICalyet et al.i (j2004f ): iFurlan et al.1 
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TABLE 3 

Properties of the targets in the Chamaeleon I region 



Name 


Spectral Type" 


Ay" 


EW(WiJ.m) 


^13-31 


Median 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


2M J10580597-7711501 


M5.25 


1.6 


3.22 


0.23 




2M J11062942- 7724586 


M6 


20.0 


1.53 


-0.27 




2M J11065939-7530559 


M5.25 


0.6 


1.70 


0.37 




2M J11070369- 7724307 


M7.5 


16.5 


2.60 


0.19 




2M J11241186-7630425 


M5 


2.7 


4.52 


0.68 




B43 


M3.25 


8.0 


6.98 


0.42 




Cl-6 


Ml. 25 


11.6 


2.34 


-0.55 


X 


C7-1 


M5 


5.0 


1.08 


-0.54 




Cha Hal 


M7.75 


0.0 


1.87 


0.00 




Cha Ha2 


M5.25 


3.8 


1.08 


-1.00 




CHSM 7869 


M6 


1.6 


1.74 


-0.57 




CHSM 10862 


M5.75 


1.6 


0.17 


-0.06 




CHXR 20 


K6 


3.5 


4.28 


-0.99 


X 


CHXR 22E 


M3.5 


1.4 


0.22 


0.63 




CHXR 30A 


MO 


10.6 


2.66 


-0.60 


X 


CHXR SOB 


Ml. 25 


11.2 


2.99 


-0.40 


X 


CHXR 47 


K3 


5.1 


2.23 


-0.66 




CR Cha 


K2 


1.5 


4.66 


-0.14 




CS Cha 


K6 


0.3 


3.02 


2.89 


X 


CT Cha 


K5 


1.6 


2.14 


-0.31 


X 


CU Cha 


B9.5 


1.5 


-0.49 


1.04 




CV Cha 


G9 


1.5 


5.78 


-0.27 




DI Cha 


G2 


2.7 


1.91 


-0.67 




Hn 5 


M4.5 


1.1 


2.07 


-1.45 




Hn lOE 


M3.25 


3.6 


3.23 


0.35 




Hn 11 


MO 


7.6 


2.56 


-0.73 


X 


Hn 13 


M5.75 


0.8 


1.29 


-0.44 




Hn 21W 


M4 


2.6 


1.45 


-0.77 




ISO 52 


M4 


1.3 


2.20 


0.11 




ISO 79 


M5.25 


7.9 


1.58 


-0.14 




ISO 91 


M3 


14.5 


7.53 


-0.40 




ISO 138 


M7 


0.0 


1.81 


-0.27 




ISO 143 


M5 


3.0 


1.47 


-1.01 




ISO 220 


M5.75 


6.1 


1.44 


-0.82 




ISO 225 


Ml. 75 


4.4 


1.11 


0.15 


X 


ISO 235 


M5.5 


7.5 


1.48 


-1.26 




ISO 237 


K5.5 


6.8 


3.43 


0.51 


X 


ISO 252 


M6 


3.4 


1.35 


-0.85 




ISO 256 


M4.5 


9.1 


3.27 


-0.63 




ISO 282 


M4.75 


3.6 


1.86 


-0.53 




SX Cha 


MO 


2.8 


2.52 


-0.48 


X 


SZ Cha 


KO 


1.5 


3.38 


1.68 




T5 


M3.25 


1.2 


1.36 


-0.32 




T21 


G5 


3.3 


-0.04 


-1.05 




T25 


M2.5 


1.6 


3.30 


2.79 


X 


T28 


MO 


4.8 


2.05 


-0.60 


X 


T29 


K6 


7.3 


2.03 


-0.34 


X 


T33A 


G7 


3.0 


5.38 


-0.54 




T35 


MO 


3.5 


0.51 


1.49 


X 


T42 


K5 


8.2 


1.46 


-0.02 


X 


T43 


M2 


5.2 


1.90 


-0.41 


X 
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TABLE 3 — Continued 



Name 


Spectral Type^ 




EW {10 iJ.m) 


Ti\^ 31 


Median 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


T45a 


Ml 


1.9 


1.66 


0.10 


X 


T47 


M2 


4.2 


2.14 


-0.28 


X 


T50 


M5 


0.4 


2.04 


-0.32 




T51 


K3.5 


0.8 


4.18 


-1.53 




T54 


G8 


1.8 


0.62 


1.10 




T56 


M0.5 


0.6 


7.57 


0.90 


X 


TW Cha 


MO 


1.2 


5.57 


-0.18 


X 


UY Cha 


M4.25 


0.0 


3.10 


-0.93 




UZ Cha 


M0.5 


2.2 


2.13 


-0.22 


X 


VV Cha 


M3 


0.5 


2.43 


-0.79 




VW Cha 


MO 


2.6 


2.52 


-0.21 


X 


VY Cha 


M0.5 


3.2 


2.58 


-0.22 


X 


VZ Cha 


K6 


2.0 


1.35 


-0.91 


X 


WW Cha 


K5 


4.8 


4.16 


0.42 


X 


WX Cha 


Ml. 25 


2.0 


2.56 


-1.03 


X 


WY Cha 


MO 


3.0 


1.61 


-1.18 


X 


WZ Cha 


M3.75 


0.5 


1.97 


-0.75 




XX Cha 


M2 


1.2 


1.77 


-0.39 


X 



Note. — Se e the footnote of Table [T] for an explanation of the table content. Uncertainties in Av, assumed to be 0.46 mag 

l|Luhmanll20"07l '). account for ~ 10% and 2% of the error bars of EW{10^Ta) and 7113- 31, resp ecti vely. 

" Spectral type and, where available, extinction information were taken from lLuhmanI l|2004al ) and lGtSroez fc MardonejipOOa) . 
For objects without published or with uncertain extinction values, the optical extinction Ay was derived using the observed 
optical/near-infrared colors and intrinsic colors based on the spectral type or a typical CTTS excess (see text for details). 



TABLE 4 

Properties of the targets in the Ophiuchus off-core region 



Name 


Spectral Type" 




EW{10iJ.m) 


"13-31 


(1) 


(2) 


(3) 


(4) 


(5) 


16126-2235 


M3 


0.6 


10.98 


1.15 


16156-2358 


FO 


2.2 


-0.02 


1.36 


16193-2314 


G5 


3.5 


3.47 


-0.14 


16201-2410 


GO 


8.0 


6.40 


1.32 


16225-2607 


K7 


1.5 


6.89 


0.52 


16289-2457 


G5 


9.5 


5.16 


0.21 


16293-2424 


GO 


17.1 


2.49 


-0.03 


DoAr 16 


K6 


3.8 


2.30 


-0.03 


DoAr 28 


K5 


2.7 


1.85 


2.19 


IRS 60 


K2 


9.0 


3.24 


-0.49 


L 1689SN02 


M3 


15.0 


0.45 


-0.83 


ROX 42C 


K6 


1.9 


2.08 


-0.37 


ROX 43A1 


GO 


3.5 


5.86 


-0.52 


ROX 44 


KO 


3.8 


8.05 


0.63 


ROX 47A 


M3 


1.9 


1.20 


-0.66 



Note. — See the footnote of Table[T]for an explanation of the table content. Uncertainties in Ay, which are typically 0.3 
mag (McClure et al., in preparation ), account for ^ 10% an d 1% of the error b a rs of i?iy(10/in i) and ni 3-3i, respect i vely. 
" Spectral t ypes were taken from iHerbig fc BeUl dlQSSl ): iChen et all (Il995h : iPreibisch etall (199^ ): iVieha et al.l (pOOl ): 
[Torres et al.l ((2006;); McClure et al. (in preparation). 
" Av values are from McClure et al. (in preparation). 



